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SYNOPSIS 


•The Mahanadi is one of the major rivers in Peninsular 

India draining into the Bay of Bengal"! It has a catchment area of 
2 

148,589 km and is occupied by a variety of rocks, alluvial 

deposits and soils. -As one of the main objectives of this thesis 

was to study the details of weathering of bed rocks, field work 

was confined to the upper part of the river basin which is 

relatively free of recent deposits and human interference. 

•The source area of the Mahanadi is occupied by granites 

and granite gneisses of Archaean age\ The river then flows for a 

considerable length through the well known Chhatt i sgarh Basin 

containing Proterozoic dolomitic limestones and calcareous shales. 

•The granitic and limestone terrains together cover 46,700 

2 

km <51.4^ of the entire drainage basin) and form the study area 
for this work*. A mi neral og i ca 1 study of rocks and soils collected 
from this area was combined with a 5-year water and sediment 
database obtained from the Central Water Commission <CWC>, 
Government of India and was supplemented by spot checks. The CWC 
sampling station# are spread over 15 locations throughout the 


Mahanadi Basin. 



Geomorpho 1 ogi cal data for longitudinal profile and 
hypsometric analysis indicate that the river is actively eroding 
granitic rocks in the source area.- Satellite imageries at the 
same time show typical karst features in the limestone terrain. 
Presence of unweathered felspar and occasional calcite in 
suspended sediments point to a transport dominated (weathering 
limited) denudation regime. The overall ratio of physical to 
chemical erosion rate works out to be around 7:1.- 

Ca^ + , Mg^ + , Na + , K + are the major cations and HCO^ , Cl , 

2 - 

50^ are the major anions in river water. The concentration of 
major ions at the confluence of tributaries matches fairly well 
with the calculated values obtained by weighting the concentration 
in individual tributaries according to their drainage areas. This 
indicates a major contribution from weathering of rocks in the 
drainage area. 

-River and ground water containing dissolved CO^ are the 

principal weathering agents: The apparent partial pressure of 

CO^- calculated from pH and HCO^ in river water, averages around 
-2 46 

10 ■ atmosphere. The total HCO, can be distributed between two 

sources namely, carbonate weathering and silicate weathering. A 
rough calculation using global averages for this fractionation 
shows that the ratio of carbonate derived HCO, to silicate derived 
HCQ^ reaches a value of 2.0 in the limestone terrain and is more 



than 1.0 all along the river course. This suggests that the 

influence of weathering of carbonate minerals is dominant over 
weathering of silicate minerals even in the downstream segments 
outside the limestone terrain. The same conclusion is supported 
by a balance between (Ca + Mg). , . and (HCO,). . . and a cluster 

of river water analyses near the HCO^ apex in a ternary plot of 

HCO, , SiO_ and SO, on an equivalent basis. An R-mode Factor 

>24 

Analysis of water composition also indicated that a Na-Ca-Mg-HCO^ 
factor can explain 70?fc of the total variance among the ions. 

'Mineral water equilibria were further verified through 
thermodynamic data on (1) carbonate and (2) silicate systems. The 
river is generally undersaturated with respect to calcite and 

dolomite which explains its aggressive weathering of dolomitic 
limestone; In the standard silicate mineral stability diagrams 
for K. Na and Mg systems. river water analyses at different 
sampling stations plot within the kaolinite field. On the other 
hand, the analyses plot near the kaolinite-Ca montmor i 1 1 on i te 
boundary in the stability diagram for the Ca-system suggesting 
that these two clay minerals can co-exist in the river water 

environment. This conclusion is consistent with the presence of 

o 

kaolinite and a glycol expansive 1 4A clay mineral in suspended 


sediments and soils. 



A limited number of analysis of heavy metals in soils 
and sediments showed a general distribution pattern of Pb < Cu < 
Mn < Zn in the granitic terrain and Pb < Zn < Cu < Mn in the 
carbonate terrain. 

A detailed study of soil profiles in a pilot area within 
the Chhattisgarh Basin confirmed the dual weathering scheme of 
carbonate and silicate minerals. While a congruent weathering of 
calcite and dolomite has produced typical solution features, 
weathering of glauconite present in s t roma to 1 i t i c limestones has 
led to the formation of a lateritic soil cover. Relative mobility 
of elements during weathering has been in the order Ca > Mg > Na > 
K > Si > A1 > Fe. 
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CHAPTER 1 


INTRODUCTION AND OBJECTIVE 

Rock weathering is a fundamental step within the fluvial 
trosion cycle. River water can cause physical breakdown of rocks 
n the drainage basin as well as chemical decomposition of the 
:onstituent minerals. As a result, the river acquires its 
sediment load and its dissolved chemical content. Therefore, a 
study of weathering processes in a river basin is important in any 
*ater quality assessment programme. 


1.1 Weathering in River Basins 

A large river basin contains enough variation in 
lithology and geomorphology so that a comparative study of 
physical and chemical processes is usually possible. The 
compilation of data on erosion rates from individual basins leads 
to a global estimate of river input to the ocean systems. Many 
fundamental geochemical concepts dealing with 


atmosphere-hydrosphere- 1 i thosphere interactions are based on such 
studies <e.g. Garrels and Mackenzie, 1971; Holland, 1978). 

A detailed investigation of a river basin may 
incorporate the following components. 

1 > Material transport 

a) Bed rock geology 

b) Geomorphology and soil characteristics 
c> Mineralogy and grain size of sediments 



d) Sediment delivery 

e> Erosion rates 
> Hydrolog i cal parameters 

a> Precipitation and infiltration 

b) Ground water recharge rate 

c> Discharge, e 1 evat i on , catchment area inter-relationships 
1 ) Geochemic a l models 

a> Mass balance of weathering reactions 

b) Thermodynamic stability of minerals 

c) Statistical evaluation 

d> Kinetics of weathering reaction 

e> Isotope studies. 


On a continental scale. the rate of physical 
(mechanical) erosion shows a wide variation. Some workers like 
Garrels and Mackenzie (1971) and Holland (1978) have observed an 
exponential relationship between this rate and the relief of 
drainage area. Similar trends are noticed among individual river 
basins for example, in Himalayan and Peninsular Indian rivers as 
discussed later. 

The rate of chemical erosion, on the other hand, varies 
within narrow limits. This rate is mainly controlled by rock 
types in the drainage basin which are subjected to chemical 
weathering under a more or less uniform C0 2 -pressure in the 


earth’s atmosphere. 



•2 Controls on River Water Composition 

Many attempts have been made to infer the source of 
dissolved components in river water. A common approach is to 
estimate the composition of primary and secondary minerals in a 
selected weathering scheme that will be consistent with the 
observed water chemistry. A simplified mass balance equation that 
depicts the input and output concepts in a river basin can be 
written as 

Precipitation + Weathering = Run-off + Storage + Biomass. 

Good examples of this line of investigation are available from 
Cleaves et al . (1970), Holland (1978). Stallard and Edmond (1983), 
Meybeck (1987), Berner and Berner (1987) and Drever (1988). 

In a study to identify the environmental factors that 
influence river water chemistry. Gibbs (1970) suggested the 
following classification of rivers. (1) Low TDS , high Na/(Na+Ca) 
waters controlled by rainfall, (2) Intermediate TDS, low 
Na/(Na+Ca) waters controlled by rock weathering and (3) High TDS, 
high Na/(Na+Ca) waters controlled by evaporation. Subsequent 
revision by Stallard and Edmond (1983) as discussed by Berner and 
Berner (1987) showed that these three groups are in fact 
representative of three types of rock weathering. For example, 
the first group is controlled by unreactive drainage basins 
containing cation-poor siliceous sediments and soils which have 
already undergone intense weathering. The second group includes 
most of the major rivers of the world and represents weathering of 
sedimentary rocks dominated by a CaCO^ component. In these 



2 + 

rivers. Ca and HCO^ would be the main dissolved ions. The third 
group, in contrast, represents weathering of evaporite beds. The 
TDS of these rivers is typically more than 250 mg/1 with 
relatively high concentrations of Na + and Cl~. In addition to the 
importance of rock types in controlling river water composition, 
Berner and Berner <1987) as well as Drever <1988) list some other 
natural environmental factors. These are 1) amount and nature of 
rainfall and evaporation, 2) temperature, 5) relief, 4) vegetation 
and 5) time. It is also obvious that in populated areas. these 
natural factors are superimposed by anthropogenic contributions 
like agricultural, domestic and industrial effluents. Meade 
<1969) has discussed the various sources of suspended and 
dissolved constituents due to human activity. 

1.3 Geochemical Studies on Weathering in Indian River Basins 

It is well known that the sub-continent of India is 
drained by some of the largest rivers of the world. Among these, 
the Ganges and the Brahmaputra jointly deliver possibly the 
largest sediment load to the world ocean <Milliman and Meade, 
1983). One of the early attempts to quantify erosion rates in the 
Himalayan region was by Raymahashay <1973>. This was followed by 
a large volume of work on Himalayan as well as Peninsular rivers 
by V. Subramanian and co-workers which has been summarized by 
Subramanian <1979, 1987) and Subramanian et al. <1987a, 1987b). 
Table 1.1 shows data from Subramanian <1979) on sediment yield and 
the apparent control of basin relief on Himalayan versus 
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Table 1.1 Sediment yield and relief of Indian river basins 
C modified from Subramanian, 1 9793 


Bas i n 

Mean elevation <m> 

Sediment yield 
(t/km^/yr) 

Himalayan Rivers 

Brahmaputra 

5000 

865 . 5 

Ganges 

3000 

591 . 5 

Peninsular rivers 

Narmada 

760 

58.7 

Tapt i 

740 

41 .7 

Cauvery 

630 

8 . 1 

Mahanad i 

500 

15.6 

Kr i shna 

420 

42 . 2 


Godavar i 


400 


56.0 



eninsular rivers. The clay mineralogy of bed loads has been 
lescribed by Naidu et al . <1985). A weathering model based on 
licarbonate content was suggested by Raymahashay (1986). More 
ecently, the major ion chemistry and the abundance of 
•adionucl ides in the Ganges-Brahmaputra system has been discussed 
>y Sarin at al . (1989. 1990). 

The Mahanadi basin specifically has been studied by 
Ihakrapani and Subramanian <1990a, 1990b). This is one of the 
Largest rivers in Peninsular India. The drainage basin contains a 
slide variety of rocks and there is some controversy about the 
lature and source of sediment mineralogy. Ray et al. (1984) had 
sarlier studied the chemical flux at the estuary of the Mahanadi 
river . 

In addition to detailed geochemical studies mentioned 
ibove, routine monitoring of all the major river basins is 
jndertaken by various government agencies like the Central Water 
Commission, Central Soil and Materials Research Station, National 
Environmental Engineering Research Institute. However, there are 
large differences in the details of water quality, sediment 
mineralogy and basin geology. There is an obvious need for a 
systematic study of the nature of bed rocks, weathering and soil 
formation with special reference to the acquisition of dissolved 
and suspended load by rivers of India. It is interesting to note 
that the Department of Science and Technology, Government of India 
has recognized some of these aspects as challenging areas in Earth 
Sciences in the next decade (Tandon and Gupta, 1990). 



1.4 Objectives of Present Work 

With this background of the present status of river 
basin studies in India, it was decided to select a representative 
area of bed rock exposures in the upper reaches of the Mahanadi 
river and go into the details of rock-water interaction. The main 
objectives of this project can be summarized as follows. 

1. To compile a mineralogi cal description of the major rock 
types exposed within the Proterozoic Chhattisgarh Basin and its 
southern fringe of Archaean basement where the source of the 
Mahanadi river is located. 

2. To distinguish between primary and secondary minerals in 
weathered rock and soil profile with the ultimate aim of setting 
up a series of weathering reactions. 

3. To establish a reliable database for the variation of water 
composition along the river course by a critical assessment of 
available literature and through spot checks. 

4. To summarize the quantitative geomorphology of the Mahanadi 
basin through a comparative study of physical and chemical erosion 
processes . 


3. To trace 

the mob i 1 i ty of 

major 

e 1 ements 

and a 

few 

environmental ly 

significant heavy 

metals 

during 

weather i ng 

and 


soil formation. 

6, To explore the relationship between rock weathering and river 
water composition by testing existing models which are based on 
Mass Balance, Thermodynamics and Statistics. 



CHAPTER 2 


SAMPLING AND ANALYSIS 


2.1 Field Work 

Field work for sample collection was confined to the 
upper part of the Mahanadi basin. Sampling of rock. soil and 
sediment was carried out from all the geologic formations exposed 
in the area. At least 3 representative fresh and weathered rocks 
were sampled from each formation. Residual soil samples were 
collected wherever a soil profile was available. Flood plain 
sediments were collected from 0-10 cm depth. Suspended sediment 
was collected at the site by filtering large volumes of turbid 
river water in the rainy season through Whatman 42 filter paper. 
Water samples were collected during post-monsoon period of 1989 
from four selected CWC stations at Baronda, Rajim, Simga and 
Andh i yarkore . A sample of monsoonal rain (August 1990) was 
collected by a specially erected rain gauge at I.G. University of 
Agriculture, Raipur. To avoid aerosols, rainwater was collected 
ten minutes after the rain commenced. pH and bicarbonate were 
determined at the sampling site while conductivity was measured in 
a nearby chemical laboratory within six hours of collection. 



2. 2 Laboratory Work 

Laboratory work consisted of chemical analysis of a 
limited number of water samples for major ions, chemical analysis 
of rocks for major oxides, mineralogi cal analysis of rock. soil 
and sediment, grain size analysis of suspended sediments and heavy 
metal analysis of sediments and soils. Prior to chemical analysis 
of water, it was filtered through 0.45 pm millipore paper and 
stored in polyethylene container in a refrigerator. Analyses were 
carried out within 15 days of sample collection. The procedures 
adopted for water analysis are summarized in Table 2.1. Bulk 
chemical analysis of a few representative rock samples were 
carried out by X-ray fluorescence (XRF) unit at Regional Research 
Laboratory, Bhubaneshwar. USGS rock standard G-l was also 
analyzed for calibration. 

2.2.1 Mineralogi cal Analysis 

Mi neral og i cal analysis of rock, soil and sediment was 
carried out by petrographic examination of thin sections. X-ray 
diffraction (Siefert XRD unit with CrKct radiation), differential 
thermal analysis (Lynseis Type 2045) and scanning electron 
microscopy (JEOL-840A model coupled with a KEVEX energy dispersive 
X-ray analytical system). X-ray studies were carried out with two 
sets of samples namely bulk powder and smear slide. Preheated 



XU 


able 2.1 Analytical methods for water samples 

pH : pH meter in laboratory within six hours of sample collection 

Bicarbonate ; Acid (HC1) titration at site 

Conductivity : Conductivity meter within six hours of sample collection 

Chloride : AgNO^ titration 

Sulphate : BaCl^ titration 

Na and K ; Flame AAS 

Ca and Mg s EDTA titration 
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powder samples were also subjected to X-ray diffraction in some 
cases. Smear slides were glycolated wherever necessary. 

2.2.2 Grain Size Analysis 

Size analysis of suspended sediments was performed by a 
Coulter Counter. Sediment suspension was prepared in NaCl 

solution to produce the desired voltage pulse which was scaled and 
counted. 

2.2.3 Heavy Metal Analysis 

Heavy metal analysis of soils and flood plain deposits 
was carried out by digestion in HNO^-HCl-HF mixture in a teflon 
bomb. The extract was analyzed in a Perkin Elmer model 2380 AAS 
unit. USGS rock standards MAG-1 and SGR-1 were used for 


cal ibrat ion. 



CHAPTER 3 


PHYSIOGRAPHY AND GEOLOGY OF THE MAHAN ADI RIVER BASIN 


The Mahanadi is the second major river in Peninsular 
India with respect to its water potential. The river basin is 
located in the states of Madhya Pradesh and Orissa. With a 
catchment area of 148589 square kilometers, the basin lies between 
the longitudes 80°30 E, 84°50 E and latitudes 19°20 N and 23°35 N. 
The Mahanadi river rises from a pool 6 kms away from Pharsiya 


village near Nagri town of M.P. at an 

e 1 evat i on 

of 

457 

metres . 

Hydrologi cal ly . the river basin can be 

d i v i ded 

into 

1 1 

smal 1 er 


sub-basins (Fig. 3.1). The principal tributaries are Pairi, 
Seonath, Jonk. Hasdo , Mand, lb and Tel. Fig. 3.2 shows the course 
of main Mahanadi and its tributaries. 

3.1 Physiographic Divisions of the Basin 

The following physiographic zones can be identified in 

the basin area. 

a) Highlands at the south-west margin. 
b> Erosional plains at western and central parts. 
c> Plateau areas towards north 

d) Eastern ghats towards south 

e) Alluvial plains towards east. 




BENGAL 








tributaries 



The highlands at the SW margin comprise parts of Durg, 
Raipur and Bastar districts in M.P. Erosional plains at western 
and central parts cover the well known Chhattisgarh region 
encompassing Rajnandgaon, Durg, Raipur, Bilaspur and a part of 
Raigarh district. Plateau areas in the northern part include 
Korba, Champa and Raigarh districts of M.P. The eastern ghats 
towards south have a gentle slope from north to south and include 
Mayurbhanj, Keonjhar, Bolangir, Phulbani and parts of Dhenkanal 
and Koraput districts in Orissa. This region is characterized by 
undulating country sides with hill ranges of 1000 metres and 
above. The coastal plains form an extensive and wide alluvial 
tract which cover Cuttack and parts of Puri district in Orissa. 
This area also contains high level laterite covered plateaus. They 
represent the uplifted and dissected erosional surfaces. These 
surfaces truncate against the eastern ghats. The mean basin 
elevation is around 500 metres. 

3.2 Climate 

3.2.1 Rainfall and Evaporation 

The overall climate is tropical to subtropical. Winter 
season normally sets over the basin from December to February. 
Precipitation during this period is mostly confined to the 
northern half of the catchment. The summer season is relatively 



long, continuing from March till June. South-west monsoon 
normally advances over the basin by the third week of June and the 
*>hole catchment receives moderate rainfall till the first week of 
October . 

There are about 212 rain gauge stations maintained by 
IMD (Indian Meteorological Department) which are fairly 
distributed all over the basin. Rainfall intensity is not marked 
by any significant variation in the Mahanadi catchment. Bulk of 
rainfall (91%) occurs in the monsoon period from mid June to first 
week of December. Average annual rainfall is 144 cm. Monthly 
distribution of annual rainfall for the entire basin is given in 
Table 5 . 1 . 

The basin has a high rate of potential 
evapo- transp i rat i on . IMD have reported that the annual potential 
evaporation varies from 174.5 cm in the west to 152 cm in the 
east. The highest relative humidity of around 80 percent is 
experienced in August. 

3.2.2 Temperature 

December is the coldest and May is the hottest month of 
the year. The western part of the basin records both the extremes 
of temperature. Diurnal range of temperature in July-August is in 
the order of 5° to 6°C and during winter season, it is as high as 



Table 3.1 


Monthly distribution of normal annual 


rainfall in the Mahanadi river basin (Source: CWC) 


Month 

Ra infal 1 (mm) 

% of annual total 

January 

14.9 

1 .0 

February 

23.4 

1 .8 

March 

17.8 

1 .4 

Apr i 1 

20 . 1 

1 .4 

May 

28.9 

2.0 

June 

214.9 

14.9 

Ju 1 y 

404.9 

28.2 

August 

390.0 

27.3 

September 

227.8 

15.8 

October 

70.6 

4.9 

November 

16.8 

1 .2 

December 

4.3 

0.3 


Monsoon total 1308.1 91.0 


Annual total 1438.1 


100.0 


JL o 


i to 16 C. The highest ever recorded temperature in the basin 
3 47.2°C at five stations Raigarh, Samba 1 pur, Titlagarh, Champa 
ad Raipur during different years while the lowest value is 3 . 9°C 
t Kanker. 

.3 Physiographic Divisions vs- Geologic Formations 

A broad correlation seems to exist between the 
opography and the geologic formations in the basin. However, it 
s difficult to evaluate any lithological control on the evolution 
if the sub-basins. Each of the physiographic zones described in 
lection 3.1 appear to correspond to a particular geologic 
J ormation. A simplified geologic map is given in Fig. 3.3. It is 
seen that the highlands at the SW margin belong to hard silicate 
,'ocks like Archaean granite and granite gneiss where as the 
srosional plains at the western and central parts of the basin 
zorrespond to Chhattisgarh Supergroup formations consisting of 
easily erodible sedimentar i es like limestone, shale and sandstone. 
Plateau areas at the northern part represent outliers of Gondwana 
formations. Sandstone, grit, conglomerate and carbonaceous shale 
are the dominant rock types in this region. The eastern ghats, on 
the other hand, constitute of Precambrian khondalite and 
charnockite suite of metamorphic rocks. The coastal alluvial 
plains consisting of eroded materials from upstream regions are of 




Fig. 3-3 Geological map of the Mahanadi basin (compiled fromGSI 
maps, 19 77 ) 




ecent age. The formations in the lower reaches of the river 
onsist of low level laterites and highly unconsolidated sands 
nth little or no clay. 

1.4 Selection of Study Area 

As discussed earlier, the present study aims at 

jnderstandi ng the weathering processes in the watershed that 

influence the chemical character of the stream water. Existing 

literature shows that most of the work on Indian rivers have been 

carried out near the river mouth. This may be due to easy 

accessibility in the coastal areas. On the other hand, it is 

obvious that any effect of rock type on water chemistry will not 

be clearly discernible near the mouth due to human interferences. 

Even the physical aspects of weathering should be more clear in 

the hard rock part near the drainage divides rather than in the 

alluvium covered lower reaches. Therefore, in this work, more 

attention has been focussed at the portion hydrologi cal ly 

considered as the upstream part of the river. It has the Seonath, 

the largest tributary, the Pairi, the Hamp and the Mahanadi river 

2 

up to 3ondhra. Together, these tributary basins cover 46, 700 km 
which is nearly one-third of the total basin area. 

The selection of the upper part of the basin was also 
prompted by the general nature of weathering as evident from the 



_andsat imagery (M.S.S.. Landsat 2, band 4, Scene No. 155-045>. 
~ * 9 • 5.4 shows that the limestone area has a typical karst 
topography giving a dark tone while the marginal part with 
non-carbonate rocks has a lighter tone. The area covered in the 
imagery includes the entire Chhattisgarh basin and the basement 
rocks. A detailed discussion of the geology of this area is given 
in the following section. 


3. S Regional Geology 

The upper Mahanadi basin contains the Chhattisgarh 
basin. This basin is one of the seven Purana basins in India 
which are believed to have developed under a dominantly 
extensional tectonic regime and fed with eroded material from the 
Archaean and Early Proterozoic continental crust of the Peninsular 
shield. They are repositories of shallow marine sediments in an 
unmetamorphosed and marginally deformed state (Kale, 1991). Rocks 
of different formations are exposed in the basin. A geologic map 
of upper part of the basin is given in Fig. 5.5. The Archaean 
group of rocks underlie the sedimentary formations of the 
Chhattisgarh supergroup. There is very little information 
available on these basement rocks which occur in the source area 


of the Mahanadi river. 




» 4 L.and<sat imagery of the upper part, of Mahanadi 


i i vnr 
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Pascoe (1950) was the first to divide the sediments in 
the basin into two broad categories viz. quartzo-conglomerat i c 
Chandarpur series and overlying Raipur series of sha 1 e- 1 i mes tone 
sequence. Dutt <1962 and 1964). on the basis of detailed field 
mapping in the southern part of the basin proposed the first 
comprehensive stratigraphic succession of the Chhattisgarh basin. 
He also indicated the possibility of correlating the same with the 
sediments of adjoining Kurnool and Bastar groups. Schnitzer 
(1975) reported the cyclic nature of sedimentation within the 
Raipur series of rocks. He identified five cycles. The 
formations within each cycle have been demarcated on the basis of 
sedimentary structures like ripple marks and stromatolites as also 
the Ca and Mg content of limestones. Murti (1987) established the 
stratigraphy of the central part of the basin. Das et al..<1989> 
and Mukherjee (1990) gave the stratigraphic sequence of north- 
western and south-western part of the basin. A comparative scheme 
of different parts of the Chhattisgarh basin adopted from 
Mukherjee (1990) is given in Table >.2. A simplified and 
integrated account of 1 i thostrat i graphy as applicable to the 
entire Chhattisgarh supergroup is given in Table 3.3. A summary 
of the mineralogy of rock samples collected from each formation is 
given in Table 3.4. 



25 



Table -3,2 Comparative stratigraphic scheme of Chhattisgarh Supergroup. 
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Table 3.4 


Summary of mineralogy of rock samples by XRD and 
mi croscopi c examination 


Formation Location Rock type Mineralogy 


Do tu 

Andh i yarkor e 

Reddish shale 

Q 

■f 

I 4 

Chi 

4 

F 

4 K 

4 Hem 

Kodwa 

Kodwa 

Grey limestone 

D 

+ 

Q 4 

C 

4 

F 





Ra j pura 

Sha 1 y 1 i mes tpne 

D 

4- 

Q 4 

I 

4 

F 

4 

K 4 

C 

Tarenga 

S i ngarpur 

Green shale 

Q 

4 

F 4 

D 

4 

C 

4 

Chi 

4 K 


S i n gar pur 

Black shale 

Q 

4 

F 4 

I 

4 

D 





S i mra 

Purp 1 e sha l e 

Q 

4 

F 4 

K 

4 

I 

4 

Chi 

4 D 

Chand i 

Sonad i h 

Purple limestone 

A 

4 

C 4 

Q 

4 

F 

4 

G1 

+ K + Chl 


Lakchanpur Nala 

-do- 

D 

4 

C 4 

Q 

4 

i 

4 

K 



Dhane 1 i 

-do- 

A 

4- 

Q + 

I 

4 

F 





Sonbarsa 

Grey limestone 

D 

4 

C + 

Q 

4 

I 

4 

K + 

I/Gl + F 


Dhabad i h 

-do- 

D 

4 

C 4 

Q 







N i pan i a 

-do- 

D 

4 

C 4 

Q 






Kha i ragarh 

Bi riayakpur 

Reddish sandstone 

Q 

4 

Hem 

4 

I 

4 

K 



Gunderdah i 

Khorsi Nala 

Pur p 1 e sha 1 e 

Q 

4 

C 4 

F 

4 

I 

4 

K + 

Chl + D 

Char mu r t a 

Ra j i m 

Grey limestone 

C 

4 

Q 4 

I 







Durg 

Purple limestone 

c 

4 

Q 4 

I 

4 

F 

4 

D 


Chandarpur 

Banrodh and 

Sands tone and 

Q 

4 

F 4 

G 1 







Mudpar 

qua r tz-ar en i t e 










Archaean 

Severe 1 

Granite and 

Q 

4 

F 4 

B 

4 

Px 




I oca t i ons 

gne i s s 











Q-Quar tz f F~Fe 1 spar , C-Ca lei te, D-Dolomite, A-Ankerite, B-B i o t i t e , I - 1 1 1 i t e 
K-Kaol inite, Mem-Hema t i t e , Px- Pyroxene, Chi- Chlorite, G1 -Glauconite 
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3. 5. 1 Archaean Basement 

Rocks of Archaean age which form the basement of the 
Chhattisgarh supergroup are exposed in the southernmost part of 
the region and covers the source area of the Mahanadi river. Grey 
as well as pink granites and their metamorphic equivalents 
constitute the predominant rock types. Grey granites are more 
abundant than the pink ones and show exfoliated surfaces. Grain 
size increases from grey to pink granites. Laxmi Narayan and 
Kuity (1988) have suggested a metasomatic origin for these 
granitic rocks. 

Thin section study shows that the granites have two 
types of felspar. plagioclase and orthoclase in addition to 
quartz, biotite and pyroxene among the major minerals. 
Plagioclase grains have pitted surfaces and rarely exhibit 
lamellar twinning indicating incipient to moderate weathering to 
clay minerals. The plagioclase composition was estimated to be in 
the range An 6 0~ An 7Q by measuring extinction angles on sections 
parallel to (001) and perpendicular to (001) and (010) 
crystallographic planes (after Wahlstrom, 1955). XRD pattern of 
the samples show characteristic peaks for an anorthite rich 
plagioclase. K-felspar grains on the other hand. are relatively 
fresh and show typical cross-hatched twinning. 
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Granite gneisses are characterized by crudely banded 
texture. They contain quartz, plagioclase felspar, pyroxene and 
iron oxides. This plagioclase is also anorthite rich and 


extens l ve 1 y 

weathered . 

Pyroxene grains 

show alteration 

to 

chlorite and 

clay. A sketch 

of typical granite gneiss in 

th i n 

section is g 

i ven in Fig. 

3.6. 

Bulk chemical 

analysis of one 

grey 

granite by 

XRF gave 

the 

following composition in weight 

percentage . 






S ‘°2 * 

69.93 


MgO = 0.38 

7S 

IN) 

O 

il 

> 

87 

A 1 ^ 0 j x 

14.81 


CaO = 1.66 



Fe 2 0, - 

4.36 


Na 2 0 = 1.66 




3. 5. 2 Chandarpur Group 

Following the nomenclature of Murti <1987) the 
Chandarpur group can be sub-divided into three formations namely 
1) Lohardih, 2) Chapordih and 3> Kansapathar/Kondkera . This group 
consists of various types of sandstone and quartz-aren i te . These 
are well developed along the southern margin of the Chhattisgarh 
basin. It is found that detrital quartz constitutes more than 90% 
of the sandstones and quartz-aren i tes . Felspars are scarce to 
absent. Occurrence of authigenic glauconite pellets in the 
sandstones of Chapordih formation and detrital glauconites in 



30 


CR 13 (75x) 

Fig. 3-6 Granite gneiss from chinnar village 
near Kanker 
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Kansapathar sandstones have been reported by Murti (1987) and Khan 
and Mukherjee (1990). 

3.5.3 Raipur group 

This group of sediments overlie the Chandarpur group of 
rocks with an erosional d i sconf ormi ty . It essentially consists of 
limestone-shale sequences with occasional surface and sub-surface 
occurrences of gypsum. The Raipur group of rocks can be divided 
into seven formations, viz. 1) Dotu . 2) Kodwa. 3) Tarenga, 4) 
Chandi . 5) Khairagarh, 6) Gunderdehi and 7) Charmuria. Some of 
the above formations are either absent or poorly represented in 
many parts of the basin. 

3. 5. 3. 1 Charmuria Formation 

This formation is well developed over large areas in the 
southern part of the basin. It is represented by grey and purple 
limestones of micritic texture. Towards Rajim, grey varieties of 
thickly bedded limestone (Fig. 3.7) are exposed while in Durg area 
massive varieties of purple limestone are noticed. 
Hi neral ogi cal 1 y Charmuria limestones are rich in calcite 
(specially in the grey varieties) and dolomite followed by 
significant amounts of quartz, plagioclase felspar, i 1 1 i te and 
kaolinite. Mukherjee and Khan (1989) recorded the presence of 
quartz, occasional mica, detrital glauconite and pyrite nodules in 
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thinly bedded 1 imestones from Lohara (Durg district). Murti 
(1987) reported a clay horizon (Sirpur member) underlying the 
Charmuria limestone in the central part of the basin. 

Dutt (1962) was perhaps the first to present the 
chemical composition of charmuria limestone (Table 3.5a) from a 
number of localities spread over a large area in southern part of 
Chhattisgarh basin (Fig. 3.8). A much recent set of analysis 
(Mukherjee. 1990) from Lohara, Durg is given in Table 3.5b for 


comparison. In both sets 

of analysis, high 

CaO 

i n 1 imestone 

corresponds with low SiD^ 

value. 

To i 1 1 ustrate 

this. 

CaO. Si0 2 , 

latitude and longitude values 

(Table 3.5a) 

were 

ut i 1 i zed to 

generate three dimensional 

trend 

surfaces by 

us i ng 

a graphics 

package SURFER. It is seen 

that 

a trough in 

Si0 2 

(Fig. 3.9) 


corresponds to a crest in CaO (Fig. 3.10). 

3. 5. 3.2 Gunderdehi Formation 

This formation is well developed in south and 
southeastern part of the basin. It is represented by purple 

shales which are medium to very fine grained, thinly laminated and 
very friable. XRD of samples from Khorsi Nala near Baloda Bazar, 
show the presence of quartz, calcite, plagioclase, 
ill ite/glauconite, kaolinite. chlorite and dolomite (Table 3.4). 
Authigenic glauconite in purple shale and authigenic baryte in 
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Table 3.5a Chemical composition of Charmuria limestone 
(Dutt, 1962) 


Samp 1 e 

Si0 2 

Ai 2°3 

Fe 2°3 

CaO 

MgQ 

Lat / Long 

C5 

11.86 

2. 

00 

47.72 

1 . 86 

21°01/80° 54’ 

C9 

20.02 

3.53 

10.85 

31.16 

0.00 

20°52/81° 17 

Cll 

31 .83 

3 . 81 

9.07 

29.59 

1 .21 

20°51/81°49 

Cl 2 

34.13 

4.02 

11.56 

27.81 

0.00 

20° 5 5/81 ° 48 , 

Cl 3 

15.09 

2.01 

15.17 

38.13 

0.00 

20° 58/81° 51 ' 

Cl 4 

23.49 

3.59 

7.78 

36 . 03 

0.00 

20° 59/82° 00 

Cl 5 

14.60 

4 . 

16 

43 . 35 

1 .20 

21° 05/82° 03" 


Table 3.5b Chemical composition of Charmuria limestone, 
Lohara, Durg (Mukherjee, 1990) 


Sampl e 

sio 2 

AI 2°3 

Fe 2 0 5 

CaO 

MgO 

k 2 o 

WSC9 

20.07 

0.80 

0.62 

42.92 

0.26 

0.06 

LC4 

14.93 

1 .01 

0.69 

45.83 

0.48 

0.21 

LI 

19-23 

1.35 

0.85 

43.15 

0.44 

0.25 

L7 

21 . 36 

3.45 

1.97 

38.21 

0.98 

0.72 

L4 

13.79 

1 .73 

1 .12 

45 . 38 

0.75 

0.24 


All values of oxides in weight percent 



Hg. J-b GEOLOGY Oh 

THE SOUTHERN PART OF CHHATTISGARH, MADHYA PRADESH. 

(after Duttj962) 
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localized grey shale of Gunderdehi formation have been reported by 
Murti <1987). Sedimentary structures like ripple marks and mud 
cracks are easily noticeable. 

3. 5. 3. 3 IChairagarh Formation 

Dutt (1964) suggested that this formation is younger 
than Gunderdehi shales. It is well developed only in the southern 
part of the basin and is represented by medium grained ferruginous 
sandstone. The sandstone is deep red in colour and contains 
abundant black iron mineral grains. Under the microscope, 
irregular black iron-oxide grains are observed to surround the 
oval quartz grains. Microscopic and X-ray analysis of samples 
from Binayakpur, 15 km south of Durg, showed the presence of 
quartz and hematite in equal proportions followed by the clay 
minerals illite and kaolinite. 

3. 5. 3. 4 Chandi. Formation 

This formation conformably overlies the Gunderdehi and 
is widely exposed in the basin. It is represented by massive 
limestones of grey and purple colour. Grey varieties are more 
fine grained than the purple varieties. One prominent feature of 
these limestones is the extensive growth of stromatolites. The 
stromatolites are either LLH (laterally linked hemi spheroi ds ) or 
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SH (discrete vertically stacked hemispheroids) type. Both LLH-C 
and LLH-S modes are common while SH-V modes are noticed in larger 
blocks of limestone. Stromatolite columns are separated by 
calcareous mud and stand out prominently on weathering. 

Structural features such as bedding planes are not very 
common although joints, both open and closed. are observed in 
massive varieties. Stylolites are found at places. Dark grey 
limestones often carry thin veins of pure calcite in them. XRD 
analysis of a number of samples from Chandi formation revealed 
that these limestones are impure. They contain silicate minerals 
like quartz, albitic felspar, i 1 1 i te/glauconi te , kaolinite and 
chlorite in addition to calcite. dolomite and ankerite. Detailed 
study of the nature of weathering of this limestone is separately 
treated in Chapter-7. 

3. 5. 3. 5 Tarenga Formation 

This formation which overlies the Chandi limestone is 
well exposed in the central part of the basin. It consists of 
shales of different colour and composition. Calcareous ones are 
generally purple with interbedded chert while the more 
argillaceous varieties are green and dark grey. X-ray diffraction 
study of shale samples collected near Singarpur village, showed 
calcite and dolomite in varying proportions while the green and 
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purple shales were found to contain additional chlorite and 
kaolinite. Dark grey varieties were almost devoid of any clay 
mineral. Montmor ill on i te could not be confirmed in any of the 
samples collected during this work although Murti (1987) reported 
occurrence of pure montmor i 1 Ionite in cherty pink shales from the 
above area. 


3. 5. 3. 6 Kodwa and Dotu Formations 

These two formations conformably overlie Tarenga shales. 
Dotu formation being the youngest of the two, has been exposed in 
the north-western part of the basin (Das et al . , 1990). Kodwa 
formation consists of reddish, medium to coarse grained shaly 
dolomitic limestone and fine grained dark grey dolomitic 
limestone. X-ray mineralogy of reddish limestone from Rajpura 
near Kodwa, shows abundance of dolomite and negligible amount of 
calcite. Other silicate minerals are quartz, illite, felspar, 
kaolinite and chlorite. Grey limestones of Kodwa, on the other 
hand, are composed of significant amounts of dolomite and calcite 
besides minor proportions of quartz and felspar. 

The lithounit representing Dotu formation is essentially 
an argillaceous red shale horizon. Extensive exposures are 
observed at Bemetra and on the banks of the Hamp river at 
Andhiyarkore (Fig. 5.11). Andhiyarkore shale consists of quartz. 
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Fig* 3.11 Bedded red shale of Dot it r . . 

-u formation expose 

Andhi yar k or e 


at 
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illite, expandable chlorite, felspar, kaolinite and hematite 
minerals as identified by X-ray diffraction analysis. These 
shales are found to contain calcite crystals in vugs. Das et al . 
( 1990 ) reported the occurrence of gypsum beds mainly within the 
Dotu formation and partly in the Kodwa formation based on few 


surface occurrence and drill hole data. 



CHAPTER-4 

PHYSICAL VERSUS CHEMICAL EROSION 


Study of river erosion and sen i + +. • 

dna sea 1 mentat 1 on processes is 

important because of their env i ronmon+o i 

environmental and economic impact. The 


effect of river on watershed i 


5 evident through physical and 


chemical weathering. A large river like Mahanadi supplies 
suspended as well as dissolved inputs to the ocean. Neglecting 
human activities, the suspended load is mainly derived by the 
physical breakdown (disintegration) of rocks in the drainage area. 
The stage of erosion is reflected in river landforms like channel 
and valley features. In a similar way, the dissolved load of the 
river is controlled by chemical reaction between river water and 
rock-forming minerals. The details of thi 
discussed below. 


»ese processes are 


4.1 Stage of River Erosion 

Longitudinal profiles, hypsometric curves and 
s e d i me nt fluxes are s ome of the mo rphological info rma t ion wh i ch 
help in expressing river erosion in a quantitative way. In order 
to assess the stage of erosion in the Upper Mahana d i basin a 
fourth order tributary near the source area was selected. A 
longitudinal profile of this stream (Fig. 4.1) was plotted by 




1 i sted 


using Toposheet No. 64H of scale 1:250.000. The data are 
in Table 4.1(a). This profile shows nick points and steep slopes 
which are indicative of an active physical erosion process in the 
upland area. The downstream segment has already approached an 
equilibrium profile. 

The hypsometric curve which investigates the 
area-altitude relationship within a drainage basin also helps in 
identifying the stage of river erosion. Following the method of 
Strahler (1952), height ratios were obtained by first determining 
the total range of contours between the confluence with the main 
river and the summit point of the Turi Nadi. The height of each 
measured contour above the confluence is then determined and 
ratios to total basin height computed [Table 4.1(b)]. The ratios 
were then plotted and a smooth curve was drawn passing through the 
data points (Fig. 4.2). It is seen from this hypsometric analysis 
that 64.4% of the land mass has already been eroded and 55.6% of 
the land mass in the drainage basin remains to be eroded. 

4. 2 Denudation Regime 

Topography produced by weathering, mass movement 
and stream run-off is complex. The concept of denudation regimes 
(Carson and Kirby. 1972) is explained in terms of weathering and 
transport in relatively flat or steep terrains. Two extremes of 
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Table 4.1<a> Data for 1 ong i tud i na 1 prof i 1 e of Turi Nadi 


E 1 eva t i on 

(m) 

Stream length 

< km) 

650 

0 

600 

0.63 

600 

1 .20 

500 

2.42 

500 

5.00 

400 

11.00 

550 

14.62 


Tabl® 4.1(b) Hypsometric data for Turi Nadi 


Absolute 

contour (m) 

Relative 

he i ght (m) 

<h> 

Total 

he i ght (m) 

<H> 

Area above 

2 

contour(km ) 

<a> 

Total basin 

2 

area(km ) 

<A> 

h/H 

a/A 

800 

450 

430 

0 

141 .036 

i 

0 

700 

550 

430 

8 . 876 

141 .036 

0.77 

0.065 

600 

250 

430 

33.566 

141 .036 

0.55 

0.258 

500 

150 

450 

63.786 

141 .036 

0.55 

0.452 

400 

50 

450 

1 32 . 426 

141 .036 

0.11 

0.959 

550 

0 

450 

141 .036 

141 .036 

0 

1 






denudation regime as discussed by Stailard and Edmond (1983) are 
a) weather j_ng_l mi^ted where rate of transport is faster than that 
of generation of weathered materials, b) transport^ imi ted where 
rate of transport is less than that of generation of weathered 
products. A weathering limited denudation regime has the 
following characteristics : 1) development of thin soil profiles 
or bare rock surfaces, 2) presence of partially weathered primary 
minerals in the stream load, 3) availability of less time for soil 
and water interaction, 4) imprint of weathering of various rock 
types on the sediment load. 

Field observations in the upper part of the 
Mahanadi basin revealed wide occurrence of limestone exposures 
without any soil cover. Thin soil profiles are developed on local 
scales. Mi neral ogi cal examination of river suspended sediments 
show the presence of unweathered felspar and occasionally calcite 
besides a host of clay minerals [Table 4.2(b)]. This may be due 
to lack of sufficient time of interaction between physically 
eroded materials and river water. Size analysis of suspended 
sediments of four upland tributaries show that about 70% of the 
material is in the silt fraction [Table 4.2(a)]. These features 
are consistent with transport dominant or weathering limited 
condition. This is further validated by a higher suspended load 
compared with dissolved load carried by the river. Details of 



fable 4.2(a) Grain size distribution of suspended sed i merits (X) 


S i ze 
( /i m) 


Sample 

1 ocat i on 


Baronda 

Ra j im 

S i mga 

Andh i yarkore 

> 20.66 

-- 

__ 

— 

— 

20.66 

2.33 

-- 

— 

— 

16.40 

6 . 30 


— 

1 . 51 

13.01 

8.75 

1 .41 

2.21 

1 .69 

10.33 

5.16 

3.77 

3 .65 

2.45 

8.20 

4.90 

3.54 

2.43 

3.72 

6.51 

21 .00 

10.70 

26.19 

27.55 

5 .16 

5.51 

16.04 

15.51 

9.13 

4.10 

5 . 84 

4.58 

4.27 

8 . 01 

3.25 

10.75 

6.70 

6.18 

5.77 

2 . 58 

6.91 

5.31 

7.73 

7.06 

2.05 

5.86 

9.24 

5.06 

8.75 

1.63 

5 . 97 

13.21 

4.35 

9.17 

1 .29 

5 .31 

15.29 

11.66 

7.27 

1 .03 

2 . 66 

4.93 

6.14 

3.92 

0 .81 

1 .72 

3.43 

3.08 

2.67 

0.65 

1 .05 

1 .88 

1 . 56 

1 . 34 

'lean s i ze (^m) 6.10 

4.33 


5.13 


= absent 


Fable 4.2(b) Mineralogy of suspended sediments 

1) Baronda : Q+I+K+F 2)Rajim :Q+I+K+F 

5) S imga : G + I 4) Andhiyarkore :Q+Mo+I+K+ Chi 

3 - Quartz: I - Iliite: K - Kaoiimte; F - Felspar; Mo - Montmor i 1 1 on i 

2hl - Chlorite; C - Calcite 
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sediment loads are discussed in the following section. 

4.3 Se-diment Flux and Erosion Rate 

In the Mananadi river. CWC has 13 gauging stations 
where measurements for discharge and suspended sediment are 
carried out on a daily basis. Therefore, available information on 
basin area, discharge and suspended load for the water year 
1980-83 were compiled from CWC sources. These were the latest 
wnen this work was started. Annual average suspended load was 
calculated from daily measurements while annual average dissolved 
ioad was estimated from monthly specific conductance measurements. 
Conductance values were discharge weighted and the average thus 
obtained was multiplied by a factor of 0.65 to get the discharge 
weighted total dissolved solids <TDS> in mg/1. Annual suspended 
loads at 13 sampling stations are presented in Table 4.3. 
Suspended ioad carried by the river during monsoon and nonmonsoon 
periods for the year 1984-85 are listed in Table 4.4. Total 

physical ioad was estimated by adding a bed load value equivalent 
to 10% of the suspended load as suggested by Blatt and Middleton 
(1972) .to the suspended load value. Estimates of physical and 
chemical erosion rates are given in Table 4.5. 

It is obvious from Table 4.4 that almost all of the 
suspended load is carried during the monsoon period. Since 

* 



51 


Table 4.3 Annual suspended load at CWC stations 


Location 


Suspended load 

(million 

tonnes ) 



1980-81 

1981-82 

1982-83 

1983-84 

1984-85 

Average 

1 ) Baronda 

8.04 

0.23 

0 . 84 

0.84 

0.42 

2.07 

2 )Ra j lm 

13.60 

1 . 90 

0.74 

0 . 59 

0.51 

3.47 

3>Simga 

3.42 

1 .70 

1 .35 

2.05 

2.38 

2.18 

4) Andh i ryarkore 

0.93 

0.39 

0.45 

0.49 

0.50 

0.55 

5 ) Jondhra 

0 . 02 

2.55 

3.70 

5.25 

8.08 

3.93 

6 >Rampur 

1 .07 

0.40 

0.91 

0.35 

0.73 

0.69 

7 ) Bamn i d i h 

7.09 

2.74 

4.96 

4.58 

7.70 

5.41 

8 >Basantpur 

36.17 

15.24 

8 .53 

16.32 

17.07 

18.67 

9 ) Kurubhata 

2.47 

1 . 58 

2.58 

2.55 

3.14 

2.41 

1 0 ) Sundergarh 

2.21 

3.09 

2.63 

2.45 

5.73 

3.22 

11 )Salebhata 

2.03 

0.64 

4.02 

0.85 

1 . 51 

1 . 81 

12 >Kantamal 

9.71 

6.39 

5.87 

5.20 

7.88 

7.01 

1 3 ) T i karpara 

43.79 

19.11 

42.13 

22.23 

13.15 

28.08 


Source; Central Water Commission <CWC> 


..--NTRAL U: RARlr 

1 | T. KA\ p dR- 


4m. M®. A. 
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Table 4.4 Seasonal suspended load distribution (year 1984-85) 


Local i on 

Sediment 

1 oad 

< tonnes ) 

Sediment 

load (%> 

Monsoon 

Premonsoon 

Monsoon 

Premonsoon 

1 ) Baronda 

419816 


-- 

100 

— 

2 > Ra j l m 

506287 


-- 

100 

-- 

3 ) S i mga 

2378124 


2403 

99.9 

0.1 

4 > Andh i ryarkor e 

504211 


439 

99.9 

0.1 

5 ) Jondhra 

8082739 


1699 

99.98 

0.02 

6 > Rampur 

729426 


-- 

100 

-- 

7 ) Bamn l d l h 

7704515 


7686 

99.9 

0 . 1 

8 ) Basantpur 

17069509 


19526 

99.89 

0.11 

9 >Kurubhata 

3141892 


8610 

99.73 

0.27 

10) Sundergarh 

5732809 


4671 

99.9 

0.1 

11) Sa 1 ebhata 

1513251 


— 

100 

-- 

1 2 )Kantama 1 

7881677 


2595 

99.97 

0.03 

1 3 )T i karpara 

1 3152194 


51441 

99.61 

0.39 

Source: Central 

Water Commission 

( CWC ) ; 

' denotes 

negl i gi bl e . 




fable 4.5 Rate of erosion 
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monsoon rainfall in the Mahanadi basin constitutes bulk (86*) of 
the total, sediment yield seems to be strongly correlated with 
discharge. Figs. 4.5(a) and 4.5(b) constructed from the data in 
Table 4.5 show that sediment load fluctuates widely with time at 
most of the sampling stations. 

At Tikarpara, the Mahanadi annually carries a total 

load of 55.24 million tons. Out of this, physical load 

constitutes 50.89 million tons while chemical load is only 4.55 

million tons. However, Chakrapani and Subramaniah (1990a) 

reported a different set of values which might be due to the 

different sampling and computational technique adopted by them. 

Erosion rate estimations for each of the tributaries (Table 4.5) 

show that the ratio of physical to chemical erosion in the entire 

Mahanadi basin up to Tikarpara is 7.1 : 1. The rate of physical 

’ 2 

erosion is the highest ( 7^7 t/ km /yr) for the upland tributary 

Pairi followed by Hasdo, Mand and lb river. This is possibly due 
to the fact that the relief of these basins is much higher than 
the base level of erosion compared with other sub-basins. If we 
consider the ratio of physical to chemical erosion rates, on the 
other hand, the values are comparatively low for Seonath and Jonk 
river basins apparently because these rivers drain chemically 
reactive sedimentary rocks like limestone. 




4.4 Grain Size and Mineralogy of Sediments 
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Sediment size distribution of suspended and bed loads 
are useful for proper understanding of transport and deposition 
mechanisms. Size distribution of suspended sediments from 
stations 1, 2, 5 and 4 during the sampling period January, 1 990 is 
given in Table 4.2(a). It is seen that 5.7 ^m fraction dominates 
over other size populations. Mean size range of suspended 
sediments is 4.55 to 6.10 /jm. As discussed earlier, nearly 70% of 


the sediment is carried in 

the 

silt 

f ract i on . 

Clay 

f ract i on 

(<2*um) constitutes 

only 

16.7% 

for 

the 

Pai ri 

river 

draini ng 

dominantly siliceous 

rocks 

compared to 

26.8% 

and 

24.4% 

for the 


Seonath and Hamp respectively. These rivers drain carbonate rich 
sedimentary terrains. Mahanadi at Rajim, on the other hand, 
carries highest percent (58.7) of sediment in clay fraction. It 
drains both silicate and carbonate rocks in its catchment. Thus 
it seems grain size is possibly influenced by the rock type in the 
drainage area. 

X-ray study of suspended sediments CTable 4.2(b)3 
reveals that quartz and ill i t e are present at all the sampling 
stations. Kaolinite. montmor i 1 1 on i te and chlorite are the other 
major clay minerals. Minor to trace quantities of felspar and 
cal cite were also identified in some samples. 



Table 4.6 


Mean size of bed sediments 


Locat 1 on 


Average mean 

diameter </.im) 

1980 - 

Pre Mon 

81 

Post 

1981 - 

Pre Mon 

82 

Post 

1 ) 

Baronda 

1 151 

1216 

1 325 

1410 

1271 

999 

2) 

Ra j i m 

1785 

2621 

2678 

2555 

2666 

2130 

3) 

S i mga 

1109 

1352 

723 

793 

1468 

906 

4) 

Andh i yarkore 

7110 

3175 

2479 

5831 

1457 

1000 

5) 

Jondhra 

1579 

1282 

2018 

1536 

1869 

1545 

6) 

Rampur 

1402 

1722 

2808 

2122 

1930 

2225 

7) 

Bamn i d i h 

709 

997 

762 

1151 

694 

687 

8 > 

Basantpur 

955 

801 

921 

902 

970 

1233 

9) 

Kurubhata 

824 

992 

1058 

1099 

1032 

1330 

10> 

Sundergarh 

1400 

1479 

1942 

1674 

2027 

1885 

11 > 

Sa 1 ebha ta 

1 1 34 

1481 

1744 

12 64 

1305 

1321 

12) 

Kantama 1 

1 316 

922 

1225 

976 

917 

1475 

n> 

T i karpara 

2159 

940 

2044 

1743 

1515 

1964 


Pre : May-June samples. Mon : October-November samples 
Post : December samples; Source : CWC 
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Bed sediment data reported by CWC (Table 4.6) for the 
year 1980-82 show that average mean size is generally higher than 
1000 ij m. Lack of significant change in grain size during 
different sampling seasons indicates that discharge is not a 
controlling factor for the size of the bed sediments at all the 
stations except Andh i yarkore . Premonsoon samples of this station, 
however, record higher grain size. 



CHAPTER-5 


RIVER WATER COMPOSITION 

5. 1 Database 

Over the years a large amount of water quality data for 
Indian rivers has been generated by many organisations and 
individuals. Among them. Central Water Commission (CWC) under the 
Ministry of Water Resources, Government of India undertakes long 
term water quality studies for major river basins of India on a 
regular basis. Water samples collected on the first of every 
month are analyzed in CWC laboratories and the data are published 
for restricted circulation in the form of separate books for each 
water year starting from June of one year to the May of the next. 

For the present study, five year CWC water quality 
database was preferred over periodic short-term water analyses. 
The reasons for this preference are the following. 

a> A continuous period of five years is reasonably long for any 
water quality assessment study. Hem (1991) suggests that a more 
intensive investigation of a smaller database can help build 
useful models compared to a wider database where measurements and 
analysis become routine. 

b) Such a database is capable of deciphering long term trends or 
lack of trends in water quality. It also takes care of seasonal 
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variation of dissolved constituents and isolates probable 
causative factors affecting them. 

c) The discharge and silt load data provide scope for better 
understanding of input-output budgets because it covers both lean 
and peak flow periods. One is likely to miss these variations and 
ultimately arrive at doubtful figures in a short term sampling 
per 1 od. 

d> The thirteen sampling stations selected are uniformly 
distributed all over the basin representing all the major 
tributaries (Fig. 3 . 3 ). 

e) Samples are collected on first of every month and are 

analyzed within a week. 

f) Analysis is carried out for sixteen major and minor chemical 
parameters e.g. pH. E.C, TDS . HCO^ , Cl, SO^ . NO.^, PO^ . NH^, SiO^, 
Ca, Mg. Na, K, Fe and A1 by standard procedures. 

Keeping the above merits in view, the CWC data can be 
safely assumed to provide a representative water analysis for the 
Mahanadi river. Therefore, in the present study, monthly water 
analyses were retrieved from CWC water-year books from 1980-1985 
for the 13 sampling stations on the major tributaries and the main 
channel. This was the latest available database when this work 


was initiated. 
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5.2 Average Water Composition 

The composition of water passing a fixed sampling 
station varies with time. A series of monthly measurements for 
five years can provide mean, extreme and standard deviation which 
summarize the record of the sampling period. Table 5.1 lists the 
maximum, minimum, mean and standard deviation for concentration of 
16 chemical characteristics mentioned earlier. The mean water 
composition for each of the 13 sampling stations has been 
separately presented in Table 5.2. This was used for subsequent 
calculations and various interpretations at different stages of 
this work. 

The mean composition of Mahanadi river water shows 
certain broad features. Like most Indian rivers, the water is 
alkaline in nature. Bicarbonate is the most dominant ion followed 
by calcium, sulphate, silica and sodium. TDS is in the range of 
76 to 343 mg/litre. 

5.3 Cat i on-Ani on Balance 

Charge balance verification of cations and anions on an* 
equivalent basis provides a check for analytical precision and 
therefore, reliability of a database. For this purpose, the mean 
concentrations of dissolved species were converted to 
micro-equivalents per litre. Charge discrepancy was then 



Tabl e 5 - 1 


Summary of CUIC water analysis data 
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t 


Stat i on 


pH 

E.C 

TDS 

F • 

A i 


no 5 

P0 * 

Sl 0 2 

Na 

K 

Ca 

Mg 

HCQj 

Cl 

SQ * 


max 

8 .22 

588 

345 

4. 66 

Q. 190 

0.20 

9.50 

0.54 

60.00 

17.00 

19.70 

47.00 

9 . 20 

155.0 

57.00 

55.00 

1 IBaronda 

rrn n 

5.49 

45 

40 

a . oo 

Q. 000 

a .oo 

0.00 

a .oo 

5.00 

2.00 

1 .00 

7.10 

1 . 20 

8.0 

4.00 

0.00 


ivg 

7.23 

102 

77 

0.57 

0 . 009 

a, 02 

0.83 

0.06 

14.02 

4.83 

2.05 

12.68 

2.88 

51 . 3 

9.71 

10.04 


>d 

6.41 

73 

52 

0.83 

0.034 

0.04 

1.76 

0.12 

9.85 

2.40 

2.80 

5.51 

1.73 

20,1 

7 . 46 

12.44 


max 

8.14 

505 

211 

4.66 

Q . 190 

0.79 

9.50 

1 .20 

45.00 

24.00 

4.95 

32.10 

7.90 

121 .0 

26.00 

52,50 

2)Ra | im 

mi n 

6.50 

47 

49 

0.00 

0.000 


0.00 

0 .00 

’2.00 

2.00 

1 .05 

7. 10 

1 . 20 

27.0 

4.00 

0.00 


avg 

7.21 

1 48 

106 

0. 52 

0.009 

0.05 

0.38 

a .os 

14.28 

7.07 

1 .93 

18.35 

3.79 

74.0 

10.05 

6.98 


sd 

0.55 

57 

35 

0. 76 

0.034 

a . 1 5 

1.44 

0.18 

8.99 

3 .95 

0.75 

5 . 84 

1 ,55 

25.3 

4,65 

9.71 


3)Simga 

max 

m i n 

avg 

sd 

8 . 3 Q 

6 . 76 
7.7 5 
0,33 

519 

1 42 
328 
91 

>52 

93 

195 

48 

2.90 
0.00 
0 . 33 
0 . 45 

0.160 

0.000 

0.006 

0.023 

2.40 

0.00 

0.15 

0 . 39 

7.70 

a . oi 

1 . 67 

1 . 42 

2.04 

0 .00 
0.89 
0.28 

49.00 

2.50 

14.12 

8.54 

42 . 50 

4 . 50 
19 . 58 

8.51 

21 .00 

1 .60 
5.01 

3 .51 

45 . 70 
18.40 
31.07 

6 . 19 

26 . 00 

1 . 90 
10.41 
4.43 

199.0 

69.0 

141.7 

33.9 

34.00 

6.00 

15.19 

5 . 46 

69.00 

3.50 

21 .10 
13.56 


max 

8 . 70 

1076 

61 3 

4 . 80 

0.160 

1 .00 

4.80 

1 . 56 

72.5 0 

120 . 00 

17.50 

64 . >0 

34.10 

318.0 

>3 .00 

164.00 

4 ) Andh i - 

m i n 

7 . 09 

1 1 6 

47 

q:qo 

0.000 

0.00 

0.03 

0 .00 

2.00 

7.50 

1 .95 

4. 50 

1 .60 

55.0 

5.00 

0 . 50 

yarkorw 

avg 

7 . 76 

579 

343 

0. 74 

0.007 

Q .06 

1 . 88 

1 .00 

19.16 

43 .61 

4.14 

39 . 56 

17 . 82 

222.2 

15.28 

57.45 


sd 

0. 54 

227 

1 1 7 

0. 94 

0.025 

0.17 

1 . 34 

0.24 

12.09 

26 . 49 

2.42 

9 . 37 

7.86 

70.9 

6.76 

29.77 


max 

8 . 54 

740 

>60 

4 . 44 

0.200 

0.52 

5 . 20 

1 .20 

40.50 

65 . 50 

1 1 . 44 

64 . 30 

26.00 

220.0 

78.00 

92 .50 

5> Jondhra 

mi n 

6.95 

81 

60 
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0 . 000 

0.00 

0 . 00 

0 .00 

2.00 

2 . 00 

1 . 6 Q 

9 . 20 

1 . 9 Q 

42.0 

4.00 

3.00 


avg 

7 . 77 

361 

213 

0.76 

0.012 

0.03 

0.15 

0.11 

16.40 

22.76 

5.09 

50.85 

12.51 

143.6 

16.55 

32.34 


sd 

0 . 47 

1 60 

72 

1.20 

0 .038 

0 . 08 

1.05 

0.21 

9.52 

12 . 54 

1 .91 

9 . 99 

5 . 99 

44.7 

11.15 

17.15 


max 

8.14 

3 6 6 

264 

2.24 

0. 120 

0.20 

4.00 

0.72 

31.50 

26 . 50 

3,20 

40.60 

11.70 

164.0 

31.00 

51.00 

6)Rampur 

mi n 

6.69 

8 5 

85 

0.00 

0.000 

0.00 

0.00 

0 .00 

0.50 

4 . 50 

0.72 

10.00 

1 . 30 

42.0 

7.00 

0.00 


avg 

7.58 

2 3 8 

153 

0 . 44 

0.007 

0.01 

1.27 

0 .07 

14.20 

13.62 

1 . 88 

26.60 

6 . 64 

114.2 

11.85 

10.35 


sd 

0 . 56 

72 

36 

0.51 

Q . 022 

0.04 

1.05 

a . 16 

8.26 

5.49 

0.54 

7.09 

2.35 

32.8 

3.82 

11.38 


max 

8.15 

431 

292 

5 . 40 

0 .160 

0. 30 

10 . Q 0 

0.52 

50.00 

25.00 

9.10 

42.00 

11.90 

113.0 

47.00 

101.50 

7 ) Bamn i d l h 

m i n 

6 . 59 

51 

47 

0.00 

0.000 

0.00 

0.00 

0.00 

2.00 

2.50 

1.45 

7.10 

1 . 20 

22.0 

4.00 

0 . 0 0 


avg 

7 . 35 

173 

126 

0 . 86 

0 . 009 

0.03 

1 .52 

0.07 

14.75 

11.16 

3.03 

17.64 

4.53 

73.0 

11 . A 0 

20.52 


sd 

0 . 38 

85 

50 

1.22 

0 . 029 

0.06 

1.86 

0.13 

8.60 

6.02 

1.61 

6.27 

1 .86 

21.9 

7.53 

21,68 


max 

8.13 

442 

252 

5 . BA 

0.160 

1.64 

9.00 

0.30 

60.00 

34.00 

5.80 

36.40 

13.90 

163.0 

31.00 

57.50 

8 ) Basant - 

mm 

6 . 58 

76 

71 
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0 . 000 

0.00 

0.00 

0.00 

0 . 50 

3 . 50 

1 . 60 

10 . 20 

1 . 90 

33.0 
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a. 00 

pur 

avg 

7 . 73 

248 

162 

0.57 

0.012 

0.05 

1 . 37 

a . 07 

16.55 

15.49 

3.13 

23 . 50 

7.73 

110.2 

12.83 

18.95 


sd 

0 .41 

90 

42 

0.84 

O . Q 33 

0.22 

1 .96 

0.17 

1 1 . 57 

8.27 

1 .10 

6.02 

3.11 

34.2 

5.31 

11.02 


max 

7.87 

1 >01 

1 089 

5.50 

0, 400 

0.50 

5.20 

0.86 

52.50 

24.50 

3.45 

20.00 

MMH 

110.0 

16.00 

73.00 

9>Kuru~ 

mm 

6.46 

>9 

44 

0.00 

0.000 

0.00 

0.00 

o.ao 

3.50 

1 .50 

1.25 

5 .10 

B8H 

18.0 

4.00 

a . 00 

bhata 

avg 

7 , 35 

124 

109 

0.96 

0.011 

0.03 

a . 47 

0.14 

17.24 

5.71 

2.22 

11.49 

ED 

55.3 

8.37 

15 . >6 


sd 

0 , 37 

1 6 9 

180 

1.29 

0.056 

0.05 

0 . 79 

0.24 

9.69 

3 . 46 

0.53 

2 . 68 

HI 

16.1 

16.06 

2.41 
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8 . 22 

214 

y — 

5 . 40 

0.250 

1.45 

27.50 

d:yr 

55.00 

22.50 

3.72 

21 . 40 

11.90 

120.0 

15.00 

74.00 

1 0 ) Sunde r- 

mm 

5.35 

43 

56 

0.00 

0.000 

0.00 

0 .00 

0.00 

1 . 50 

2.00 

0.95 

7.50 

1.20 

22.0 

3. 00 

a . oo 

garh 

avg 

7 . 39 

129 

94 

Q .73 

0.008 

Q . 06 

1.21 
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15.91 

9.76 

1 .48 

13.64 

4.02 

67.7 

8.70 

13.28 


sd 

0.49 

4 3 

21 

1.19 

O . Q 35 

Q .21 

3.59 

0.18 

10.30 

4.21 

0.52 

3.27 

1.90 

20.1 

2.95 

13.11 
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8 . 40 

412 
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2.32 

0.160 

0.40 

2.70 

1.31 

28.00 
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13.20 

189.0 
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mSaia- 
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6.97 
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1 .00 

12.20 

1 .20 

41.0 

5.00 

0.00 

bhata 
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7.83 

281 

175 

0.42 
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0.05 

0 . 71 , 

0.06 

14.65 

22.32 

1.76 

26.39 

7.90 

137.5 

11.86 

9.00 


sd 

0.33 

85 

35 

0.49 

0.023 

a. 09 

0.68 

0.19 

6.56 

7.89 

0.45 

6.04 

2.70 

38.5 

4.30 

11.84 
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>60 

230 

4.76 

I 
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1 .24 
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2 ? . 00 

52 . 00 
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94 
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1 .90 
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4.00 
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rna 1 
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7.64 
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10.42 


sd 

0.43 
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416 
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1 . 8 50 
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1 . 34 
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1 . 50 
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7.47 
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0. 040 
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2. 21 

1 . 00 
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20.18 
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87.5 
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9 . 33 


sd 

0 . 36 

59 

36 

1.11 

0.239 

0.07 

2.72 

0.21 

7.74 

2.65 

0.45 

4 .43 

2.30 

32.3 

3.11 

1.21 









Table 5*2 Average water composition of the Mahanadi at different stations. 
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E.C (Electrical conductivity) in micro-mhos/cm, all other parameters except pH are in mg/ I . 
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estimated for water of each stations using the formula 


D 


£l 4 ] - [TZ ~ ] X 100 
ITl 3+CTZ - ] 


where TZ is total cations, TZ is total anions and D is percent 
of discrepancy. Table 5.3 shows that the maximum discrepancy is 
8.4 percent which is well within the tolerable limit. A plot 
(Fig. 5.1) of TZ + versus TZ resulted in a linear regressi onal 
relationship between them. The equation of regression line by 
best fit method is TZ (/ueq/1) = 1.11 TZ + ( j ueq/l) - 258.75. 

Before finally adopting the CWC data, spot checks for 
water composition were carried out by collecting samples at a few 
CWC stations located in the upper part of the Mahanadi basin 
during the field season of January-February , 1989. These waters 
were analyzed for some major ions like HCO^, Ca, Mg, Na, K and Cl 
so that an order of magnitude of reported values could be 
assessed. It was seen that most of the parameters analyzed during 
the present study fell within the concentration ranges reported by 
CWC (Figs. 5.2). 

5. 4 Expression of Water Analyses 

Various terms and units are commonly used for -reporting 
water composition. Table 5-2 shows average water analyses in 



Table >.5 Cation - anion balance 
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Fig. 5.1 Balance of total cations CT2 + 3 and total anxons CTZ 1 
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Fig. 5. a CWC vs Present water analysis (Baronda) 
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CWC vs Present water analysis (Simaa} 
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parts per million (ppm) at 13 sampling stations. It is assumed 
that ppm values are equal to milligrams per litre because dilute 
river water has a density value practically equal to 1 gm/cm^. 
Conversion of mg/litre values to mole/litre, estimation of ionic 
strengths, activity co-efficients and ion activities, were carried 
out by a user friendly computer program (Appendix I) written in 
TURBO-PASCAL using the following relations : 


mole/litre (m.) 

i 


mg/ 1 j tre x 1Q -3 

formula weight 


( 1 > 


where 'i' is any dissolved species. 


ionic strength (I) 


1 ^ 7 Z 
— Z m . • Z . 

2 l l 


( 2 > 


where Z. is the charge of the species. 

log r . = - A.Z? [j£rr - °- 2! ] °> 

where r is the activity co-efficient of the species and A is 

i 

dielectric constant equal to 0.5085 at 25 C. 



Table 5.4 Concentration of dissolved species (rtu 1 1 i -mol es/l i tre) 
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Table 5.5 


Ionic strengths and activity coefficients 
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Station Ionic Ionic activity coefficient 

strength Univalent Bivalent Trivalent 


1 ) 

Baronda 

0.0018 

2 > 

Ra j l m 

0.0023 

5) 

S i mga 

0 . 0053 

4) 

Andh i yarkore 

0.0077 

5 > 

Jondhra 

0.0053 

6) 

Rampur 

0.0036 

7) 

Bamn i d i h 

0.0028 

8) 

Basantpur 

0.0037 

9) 

Kurubhata 

0.0020 

10) 

Sundergarh 

0.0023 

11 ) 

Sa 1 ebhata 

0.0040 

12) 

Kantama 1 

0.0031 

15) 

T l karpara 

0 . 0028 


0.9539 

0 . 8278 

0.6537 

0.9483 

0.8088 

0 . 6203 

0.9248 

0.7313 

0 . 4946 

0.9115 

0.6904 

0.4345 

0.9247 

0.7313 

0.4946 

0.9366 

0.7697 

0.5549 

0 . 9434 

0.7923 

0 . 5923 

0.9359 

0.7672 

0 .5509 

0.9516 

0.8203 

0.6405 

0 . 9483 

0.8088 

0 . 6204 

0.9336 

0.7597 

0.5 388 

0 . 9408 

0.7834 

0.5774 

0.9435 

0.7924 

0.5923 



Table 5.6 Activity of dissolved species (x 10 



15) Tikarpara 0.523 
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here a t is the activity of the ion. Table 5.4 lists conversions 
n mole/litre. Ionic strengths and activity co-efficients are 
iven in Table 5.5 while ion activities are presented in Table 
. 6 . 

.5 Seasonal Variation 

Effect of season on physical and chemical parameters is 
mportant to the assessment of surface water composition. One of 
he advantages of a long term database lies in its ability to 
ecord seasonal variations, if any. Based on the distribution of 
ainfail (Table 5.1), three seasons can be identified in the 
lahanadi basin. They cover the periods 1> January to May <7.6% of 
otal rainfall), 2) June to September (86.2%) and 5) October to 
iecember (6.2%). These periods can be termed as 'pre-monsoon', 
monsoon' and 'post-monsoon' respectively. Accordingly, all major 
larameters and discharge values from the CUIC raw database have 
ieen averaged. These are given in Table 5-7. 

Certain broad features are easily discernible from 
:he diagrams (Fig. 5.3) constructed using the data in Table 5.7. 

S s usual, the effect of dilution is observed as a marked decrease 
. n the conductivity value. ’ pH is largely constant during 
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ib 1 e 5 . 7 Seasonal variation of water composition 


5 tat ion 

Season 

D i scharge 

pH 

E .C 

Si0 2 

Na 

K 

Ca 

Mg 

hco 5 

Cl 

S0 4 



Pro 

0.33 

7.44 

135.7 

11.1 

5 . 3 

2 . 0 

14.0 

5.3 

58 . 4 

11.7 

4.6 

> 

Baronda 

Mon 

61.63 

7.08 

86.7 

18.5 

4.4 

1 . 8 

12.0 

2 .0 

47 . 7 

8.4 

16.0 



P os t 

14,05 

7.13 

06.9 

11.6 

4. 5 

2 . 6 

12.0 

3 . 5 

46 . 6 

9. 1 

10.6 



Pro 

2.99 

7.13 

199.7 

15.7 

9.9 

2 . 4 

23 . 4 

4.8 


■HE9 


> 

Ra | i m 

Mon 

95 . 30 

7.09 

103.6 

17.6 

5.2 

2 . 1 

13.9 

2.8 

54 . 4 





Post 

30.61 

7 . 37 

130.8 

10.6 

5.8 

1 . 6 

16.4 

3 . 6 

69 . 1 

ill 




Pr e 

10.70 

ESS 

592.2 

12. 1 

24.4 

6 . 2 

32 . 9 

12.6 

161 .9 

18.0 

19.7 

) 

S i mga 

Mon 

144.97 

EBr 

252.0 

17.4 

16.6 

5.0 

22.9 

7.9 

111.6 

13.3 

29.7 



Post 

143.62 

HI 

321.0 

11.7 

15.3 

2 . 2 

35 . 5 

10.1 

149.0 

12.7 

10.9 



Pr a 

2 . 06 

CO 

64 2.8 

16.3 

54 . 4 

4.3 

42.9 

22.4 

247.2 

16.6 

59.6 

) 

Andh i - 

Mon 

1 1 . 38 

7.65 

480.9 

24.7 

37.4 

4.3 

34.1 

12.6 

180,8 

13.6 

59.5 


yarkore 

Post 

3 9,65 

7.85 

499.2 

15.3 

29.7 

2.9 

41 .0 

16.6 

210.8 

13.9 

48.5 



Pr e 

14.63 

7.91 

402.7 

14.2 

31.6 

4.5 

55 . 0 

16.1 

170.1 

20. 1 

30.5 

> 

Jondhr a 

Mon 

359.63 

7.34 

278. 3 

20 . 5 

18.8 

4 . 1 

26.0 

10.7 

112.4 

22.3 

32.6 



Pos t 

204. 30 

7.75 

511.9 

13.9 

14.8 

2 . 5 

34 . 7 

10.1 

145.1 

13.5 

16.2 



Pi a 

0.41 

7.76 

290 . 5 

1 1 . 2 

14.7 

1 . 8 

31.5 

8.5 

137.8 

10.7 

5.8 

) 

Rampur 

Mon 

07. 35 

7 . 36 

16 7.5 

18.2 

10.9 

2.0 

18.4 

4.9 

86 . 4 

12.9 

16.9 



Post 

2 9.70 

7 ,68 

248. 5 

11.5 

13.7 

1 . 6 

28 . 5 

6 . 6 

117.9 

11.2 

B .6 



Pro 

1 4 . 30 

7.4 3 

200 . 2 

11.9 

12.2 

3 . 3 

20.0 

5.2 

82.6 

13.7 

15.7 

) 

Ramn i d i h 

Mon 

410.60 

7.12 

157.3 

19.7 

10.7 

3 , 1 

14.1 

3 . 2 

58.7 

9.7 

40.4 



Pan t. 

76 . 06 

7,40 

151,6 

12.3 

0 . 5 

2.4 

16.6 

5.1 

71.2 

10.0 

8.6 



P r a 

55.86 

7,80 

296,3 

14.3 

20.1 

5.5 

24.9 

9.0 

123.4 

13.8 

18.2 

) 

Basant - 

Mon 

1253.60 

7.3 3 

194.1 

21 .4 

13.6 

3.2 

19.4 

6.4 

86.8 

12.0 

24.5 


pur 

Pos t 

806.20 

7.82 

2 36.0 

13.7 

10.8 

2 . 5 

26.9 

7.0 

114.0 

11.6 

12.3 



P r e 

6.31 

7 .52 

17 5.4 

14.1 

7 . 0 

2 . 4 

12.3 

4.1 

53.4 

9.3 

6.5 

> 

Kuru- 

Mon 

206 . 10 

7.02 

73.9 

22 . 6 

3.7 

2 . 2 

9.2 

3.1 

44 . 5 

7.0 

31.2 

j, , 

bha t a 

Pos t 

4 5.40 

7.4 5 

109.2 

15.5 

5 . 1 

2 . 0 

12.4 

4 . 2 

61.0 

8 . 2 

11.0 



Pr e 

5 . 50 

7 . 59 

"l 55 . F 

“1 3 7 ~ 

~12?6 ~ 

“l .4 

“I 575“" 

4 ; 4 

“7 67“ 

““10. 6 ~ 

~n 

0) 

Sunder- 

Mon 

296.20 

7.05 

92.8 

20 . 5 

7.4 

1 . 8 

11.4 

2.4 

52 . 2 

7. 3 

26.3 


gar h 

Post 

37.90 

7 , 56 

138.5 

12.7 

8 . 9 

1 . 2 

15.7 

5.8 

74 . 6 

8. 5 

14.5 



Tv#T“ 

1 .79 

7.94 

348.6 

15.5 

27.3 

1 .7 

30.4 

10.0 

167.6 

14.0 

4.7 

1) 

Sa 1 1 *« 

Mon 

15 7.00 

7.55 

179.7 

17.5 

17.7 

1 . 9 

19.3 

4.5 

90 . 8 

9.9 

20.8 


bha ta 

Post 

31.28 

7.91 

270.8 

11.9 

19.2 

1 . 6 

27.0 

7.9 

135.3 

10,1 

4.4 



m p~— 

6.00 

7.82 

24 3.2 

14.4 

17.8 

2.4 

22.1 

7.1 

111.6 

11.4 

4.4 

z > 

Kant a - 

Mon 

607 . 20 

7.27 

159.4 

24 . 0 

12.0 

2 . 2 

17.4 

4 . 8 

83.2 

9.0 

19.4 


ma l 

Pos t 

170.80 

7.75 

209.8 

17.5 

11.8 

1 . 9 

21.8 

6.5 

100.3 

11.7 

13.1 



Pr e 

~ 6780 ~ 

7.63 

201 . 5 

11.9 

9.1 

1 . 6 

18.8 

5 . 7 

"Too To” 

10.7 

5.1 

3 > 

T i kar - 

Mori 

3048.20 

7.24 

149.0 

19.8 

7.5 

2 . 1 

18.6 

4 . 2 

77.9 

9.5 

14.1 


para 

Pos t 

1445,00 

7.51 

101.4 

11.9 

6 . 5 

1 . 5 

18.9 

4 . 7 

87 . 3 

12.2 

11 .5 

i s charge in cubic 

met res/sec, E , 

C (Electric Conduc t i v i t y ) 

in micro-mhos/cm, 

a 1 1 other 


onctntrat int i n mgs/1 i tre except pH 


STATION 

O POST MONSOON 


PRE MONSOON 


A MONSOON 



1 I i riii 

2 3 4 5 6 7 8 

STATION 

PRE MONSOON O POST MONSOON 


A MONSOON 










PRE MONSOON 


STATION 

POST MONSOON 


A MONSOON 



CALCIUM(ppm) 


81 



STATION 

□ PRE MONSOON O POST MONSOON A MONSOON 




MAGNESIUM(ppm) 






BICARBONATE(ppm) 






84 


pre-monsoon and post-monsoon periods but is relatively low during 
peak flow periods. Calcium, magnesium and bicarbonate also show a 
similar trend of lowest concentrations during monsoon where as 
highest concentrations are observed in the lean season. 
Post-monsoon period is characterized by moderate values for the 
above ions. High values during monsoon period may be attributed 
to ground water recharge which is a usual feature in carbonate 
terrains. Sulphate shows highest concentrations during monsoon, 
invariably moderate during post-monsoon and low in the dry period. 
However, the Seonath river waters at the upstream region (stations 
4 and 5> do not follow this trend. Such notable deviation may be 
due to the bed-rock geology. This part of the basin consists of 
evaporite deposits at the subsurface and surface levels (Chapter 
3) . 

Silica seems to be directly related to discharge. It 
consistently maintains high concentrations during monsoon. Davis 
(1964) suggests that this may be due to rapid dissolution of 
silica during storm run-off which exposes much larger surface area 
of suspended silt and clay size particles than during the periods 
of low discharge. In addition, much run-off, during and shortly 
after the monsoon, passes through the upper part of the soil 
profile which may contain silica leached from the soil. 
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5.6. Spatial Variation 

An attempt has been made to find out whether the 
Mahanadi waters show any definite trend in major dissolved 
constituents with respect to location of sampling stations. Two 
pathways of water flow through sampling stations 1, 2. 8, 13 and 
3. 5, 8. 13 along the downstream direction of the Mahanadi river 
channel were considered. Average concentration of two groups of 
dissolved species Ca-Mg-HCO^ and Na-K-Cl-SO^-SiO„ along these 
routes have been examined for any possible spatial variation. 
Figs. 5.4 and 5.5 plot the concentration of the above species 
against distance downstream. It is observed (Fig. 5.4a,b> that 
water compositions at stations 1. 2 and 13 do not reflect any 
marked change. It is seen from the geologic map ( F i g . 3 . 3 > that 
waters at stations 1 and 2 drain mostly crystalline rocks in their 
catchments. Station 13 is located farthest downstream where 
crystalline rocks are dominant. This partly explains their 
similar water compositions. At Basantpur (station 8), all the 
constituents registered highest level of concentration. The 
increase in Ca , Mg and HCO^ values may be due to contribution of 
the Seonath tributary which drains predominantly carbonate and 
associated rocks. In addition, Jonk and Hasdo tributaries which 
drain mostly crystalline rocks, could possibly supply silicate 
derived ions and SiO^. 



Concentration (ppm) 





Spatial variation of Ca, Mg and HC0 3 along 


sampling stations 1, 2, 8 and 13 


Fig. 5. 4Ca) 
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Distance downstream (Kms) 


Spatial variation of Na, JC, Cl, SO^ and SiO^ 


along sampling stations 1, 2, 8 and 13 


Fig. 5. 4.C b> 
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Figs. 5.5(a) and 5.5(b) show that at Jondhra (station 5). 

the river water records highest concentration of Na, Cl, S0 4 and 

SiO£ wh 1 le Ca , Mg and HCQ^ concentration levels are maintained 

with respect to those at Simga (station 3). Further downstream, 

concentration of all the dissolved constituents except Si0 o 

decrease. The sharp rise in the concentrations of Na and SO at 

4 

station 5 may be due to the contributions from the Hamp tributary 
which has high concentrations of Na and SO^ compared to waters 

from other parts in the Mahanadi basin. 

5.7 Correlation Analysis 

Correlation ana lysis is a statistical tool to help 
establish some degree of association among a set of variables in a 
manner unaffected by their measurement units. In this analysis, 
correlation co-efficient (r) which is the ratio of the co-variance 
of two variables to the product of their standard deviations, is 
estimated for all the variables. This was carried out by taking 
mean values of 16 chemical parameters for all the sampling 
stations listed in Table 5.2. Computer processing of this data 
with SPSS (Statistical Package for Social Sciences) resulted in a 
correlation co-efficient matrix (Table 5.8). As seen from the 
matrix, a strong positive correlation (r = 0.995) exists between 


electrical conductivity and total dissolved solids. 


Sodium, 



Table 5,8 Correlation coefficient matrix {confidence level 33 %) 
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calcium, magnesium and bicarbonate seem to be strongly correlated 
among themselves. The correlation co-efficient is greater than 
0.9 in each pair of ions, thus indicating a common source for 
them. 


5-8 Water Quality-Discharge Relationship 

It is well established that the concentration of 
dissolved constituents in water is strongly related to river 
discharge. A cons. stent relationship is expected between average 
water discharge and solute concentration of streams in a 
conceptual dilution or mixing model. Dissolved solute 
concentration of a set of tributary waters add up at a downstream 
point. Solute concentrat i on at this point represents the diluted 
version of ail the earlier waters. If discharge factor is 
significant, then the predicted value provides a good comparison 
with the analytical value. Predicted solute concentration of 
mixed water derived by dilution can be given by the equation 


Q 2 C 2 


Q C 
n n 


n+1 


where C, , C„ C are average solute concentrations 

i Z n 


different tributaries at stations 1, 2 ....n where Qj , 
are corresponding discharge values. 


the pr ed i ct ed 
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concentration of mixed water whose discharge value is Q 

n+1 ' 

The above model was put to test in the upper part of the 
Mahanad i river. Waters of stations 2, 5. 6 and 7 mix up near 
sampling station 8 (Fig. 3. 3). Mean discharges for the above 
stations were calculated from CWC water quality database. Average 
water compositions of 2, 5, 6 and 7 were used to predict the 
composition at station 8. A bar diagram (Fig. 5.6) was 
constructed based on the results of the calculations presented in 
Table 5.9. It is seen that observed concentrations of dissolved 
species such as Ca . Mg, HCQ^, Na , K , Cl, SO^ and SiO^ is generally 
higher than the predicted values. Observed TDS is 20% higher than 
the predicted value. This surplus TDS of 20% (which is largely 
due to the bicarbonate species) may be derived from the weathering 
of carbonate rocks in the upper part of the Mahanadi catchment 
area . 

5.9 Wat*»r Quality-Catchment. Area Relationship 

A broad approach towards deciphering the relationship 
between water quality and catchment area is based on a principle 
similar to that outlined in the earlier section. Such a model 
takes into account the area weighted average water composition of 
tributaries to predict the solute composition at a downstream 
point of the main river. IT only lithologic factors are 



Table 5« O Est.imat.ion of" Q,— weighted concentration CppnO of" solutes 

in mixed water 


Dissolved F 

solute 

3 redicted concentration 
[ ° 2 C 2 + Q 5 C 5 + Q 6 C 6 + D 7 C 7 1 

[ Q 8 j 

Observed 
| concentration 

1 

Ca 

19.84 

23.50 

Mg 

6.79 

7.73 

HCCL 

> 

88.59 

110.23 

Na 

13.27 

15.49 

K 

2.39 

3.13 

Cl 

10.98 

12.83 

S0 4 

19.29 

18.95 

sio 2 

12.32 

16.55 

TDS 

173.47 

208.41 


Note 


: Q , Q , Q,, Q_ and Q„ are average discharge values at 

2 ? 6 / o 

stations 2. 5. 6, 7 and 8 respectively where C £ , C ? . C & , 

r , and C„ are corresponding solute concentrations. 

/ B 


Mg HC0 3 Na 


K Cl S0 4 Si0 2 TDS 


Predicted Cone. KsXvn Observed Cone. 


Fig. 5.6 Comparison of Q-weighted concentration with 


observed concentration 
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significant, then this predicted value provides a good comparison 


with the analytical value 


Predicted solute concentration of 


mixed water derived from upstream tributaries can be given by the 


equat i on 


A 1 C 1 + fl 2 C 2 + 


+ A C 
n n 


where C. 


. . are average solute concentrations of 


tributaries with catc hme nt areas A., A _ A respectively. C 

12 n p 

is the predicted solute concentration at the downstream station 

where the catchment area is A 

n+1 

In the upper part of the Mahanadi basin (Fig. 3. 3), 
waters from upper Mahanadi, Seonath, Jonk and Hasdo tributaries 
add up near Basantpur (station 8). Average solute concentrations 
at stations 2. 5 , 6 and 7 of above tributaries were used along 

with their catchment areas to predict the concentration at station 
8. Fig. 5.7 based on the calculations listed in Table 5.10, shows 
that observed concentration of most of the rock-derived cations 


like Ca , Mg, Na and K are equal to the predicted values. 
Magnitudes of observed bicarbonate and silica are found to be 
slightly higher while sulphate is lower than the predicted ones. 
Analytical TDS (200 ppm) is about A% higher than the area weighted 
TDS (208 ppm). Therefore, it seems that the theoretical area 
weighted water composition is a fairly good approximation of the 


analytical one. 
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T^bl e 5* 1 0 Esbx mati on of* ar ©a woi r»v*+ - j . . , . , , , 

area weighted predicted concentration 

tppnO of solutes in mixed water 


Dissolved Pi 

solute 

dieted concentration 

r a 2 c 2 4- a 5 c 5 + a 6 c 6 + a 7 c 7 ] 

L 

Observed 
| concentration 

! 

Ca 

22 . 92 

23.50 

Mg 

8.09 

7.73 

HC0_ 

3 

102.94 

110.23 

Na 

15.32 

15.49 

K 

2.90 

3.13 

Cl 

12.49 

12.83 

S0 4 

21.63 

18.95 

Si0 2 

13.78 

16.55 

TDS 

200.07 

208.41 


Note : A,, A _ , A,, A_ and A Q are catchment areas (up to stations 
2.5,6 ,7 and 8) whose values are 8760, 29645, 2920, 9730 and 
57780 sq. kms . respectively. and C Q are 

correspond i ng average concentration of solutes at stations 


2. 5. 6, 7 and 8. 
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CHAPTER-6 


ROCK-WATER INTERACTION 


6.1 Introduction 

Chemical weathering of rocks is an expression of 
rock-water interaction at relatively low temperature and pressure 
which prevail at the earth's surface environment. Extensive field 
and laboratory work has established that natural water is the 
principal weathering agent for rocks. Commonly observed 
rock-water reactions can be classified into three major groups. 
These are i > Congruent dissolution of carbonate rocks, 

1 i ) Incongruent decomposition of silicate rocks, iii) Oxidation- 
Reduction which is important in specific cases. 

The rocks exposed in the upper part of the Mahanadi 
basin offer an opportunity to study all three types of reactions 
between rock and water. Following the approaches listed by Drever 
< 1988 ), available data on mineralogy and water composition were 
used to test existing models based on Mass Balance, Thermodynamics 
and Statistical Methods. 

6.2 Role of Carbon Dioxide 

It is interesting to note that much of the weathering 
capability of natural water is due to the mild acidity produced by 
dissolved carbon dioxide. The principal source of CO^ is the 
earth's atmosphere with additional supply from decay of organic 
matter in the soil zone. The chemistry of dissolved carbon 
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dioxide can be described with the help of the following reactions. 

a) Solubility of carbon dioxide : 

co 2 + h 2 o = h 2 cd, 

K C0 2 = a H 2 C0j/ P C0 2 (1> 

b> First dissociation of carbonic acid, H^CO^ : 

* H + + HCOj 

a H + * a HC0~ / a H 2 C0 5 <2> 

c > Second dissociation of HCO^: 

HCO~ = H + + C0 2 ~ 

K 2 = a H + ' a C0 2 ”/ a HCO~ <?> 

The equilibrium constants of the above reactions have values 

* -10" 1 * 5 , K, * 10 -6 * 4 and K. = lo" 10 ' 5 at 25°C (after Carrels 

C0_ I L 

£. 

and Chr 1 st , 1965). 


H„CO_ 

t ;• 


K. 


Combining equations ( 3 ) and ( 2 > 


a, ,+ 
H 


.a 


HCO, 


10 


- 7.9 


X P 


CO, 



101 


Taking logarithm, 

109 P C0 2 = 7 -9 + log a HCQ - - pH <4> 

Equation 4 can be used to calculate the value of partial pressure 
of CO., from pH and bicarbonate content of river water at an 
average temperature of 25 C. This CO^ pressure can be defined as 
'apparent 00^, ’ of river water consistent with the given 
temperature, pH and bicarbonate alkalinity. 

The results of this calculation for the sampling points 
along the Mahanadi river have been given in Table 6.1. Log a 

HCD 

was calculated from the analysis listed in Chapter 5 . These 
calculations show that the apparent CO^ pressure in each of the 

tributaries and the main Mahanadi river is higher than the 

- 5 5 

atmospheric = 10 ' atm. with an average log =-2.464 ± 

0.1. This is in agreement with the global trend described by 


Garre 1 s and 

Mackenz i e 

<1971) 

where average 

p co 2 

of 

major world 

rivers is 

about ten 

times 

higher than 

that 

of 

the earth's 

atmosphere . 

Ho 1 1 and 

(1978) 

and Stumm and 

1 Morgan 

(1981) have 


discussed the possible reasons for the apparen 

1 Pc0 2 being higher 

than the atmospheric value. Two possibilities are 1) A perennial 
river like Mahanadi contains a significant fraction of high CO^ 
ground water, 2) The rate of release of excess C0 £ for 
re-equi 1 ibrium with the atmosphere is slower than the rate of 

solubility of CO^ gas. 

Fig, 6.1 which is a plot of apparent P CQ ^ versus TDS 
shows that the apparent CO^ pressure remains fairly constant all 
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FIG 6.1 Variation of apparent Pec* with IDS 
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over the basin, independent of dissolved load along the course of 
the river. On the other hand, subtle variations from station to 
station indicate that samples from the upper part of the basin 
where limestone and calcareous shales are dominant (Fig. 3.4) have 

relatively higher P CQ . 

6.3 Silicate versus Carbonate Rocks 

As discussed earlier, the Chhattisgarh basin occupying 
the upper part of the Mahanadi river is dominated by dolomitic 
limestone with intercalations of calcareous shale and sandstone 
(Fig. 3.4). The basin is fringed by basement rocks of granitic 
composition near the source of the river. On the other hand, in 
the lower reaches, the river passes through sandstone, granite, 
khondal l te and other metamorphic rocks which are covered by 
lateritic soils and recent alluvium. Thus it will be reasonable 
to exoect that a dual weathering scheme of carbonate and silicate 
rocks will influence the river water composition. 

6. 3. 1 Sources of Ca and Mg 

The main source of Ca and Mg in river water is the 
weathering of rock~f ormi ng minerals. Minor amounts can be derived 
from cyclic sea salts and pollution (Berner and Berner, 1987). 
The important Ca-bearing' minerals are calcite and dolomite in 
limestones; gypsum and anhydrite in evaporite beds; p lag i oc lase 
and other silicates in igneous and metamorphic rocks. Similarly, 
Mg in river water can be derived by the weathering of dolomite. 
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biolite. amphibole, pyroxene and olivine. 

Holland (1978) attempted to estimate the contribution of 
Ca and Mg from rock weathering by a comprehensive review of rock 
and river water composition. He concluded that 86% of the TDS in 
average river water is derived from sedimentary terrains. : On the 
other hand. 85 ± 8% of total Ca and 50 ± 20% of total .Mg in 

sedimentary rocks are present in carbonate minerals. Details of 
this estimation are presented in the Appendix-II. It follows that 
in r l ver wafer . 

Ca from carbonate minerals 

= <0.85 ± 0.08) X 0.86 X Ca^ . , 

total 

= 0.75 t 0.08 £ 74 ± 10% 

Mg from carbonate minerals 

= <0.50 ± 0.20) X 0.86 X Mg t j 

= 0.45 ± 0.20 ^ 40 ± 20% 

The remainder of Ca and Mg in river water is derived largely from 
silicate minerals, neglecting the minor contribution of sulphates 

and other evapor ite minerals. 

Meybeck (1987) arrived at a similar distribution of Ca 
and Mg in river water on a global scale based on a study of 
chemical weathering in selected rivers of France. His 

distribution of Ca was 67% from carbonates. 26% from silicates and 
7% from evapori tes . The values for Mg were 42% from carbonates, 
48% from silicates and 10% from evaporites. It is obvious that 
this distribution may change from basin to basin depending on the 
proportion of different rock types. However. Holland’s original 
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estimate that 74% of Ca and 40% of Mg in river water are derived 
from carbonate weathering appears to be a reasonable 
approximation. Therefore. this fractionation was adopted for 

estimating the source of bicarbonate as discussed in the next 
sect i on . 

6. 3. 2 Source of HCO^ 

The main sources of HCO ^ in river water are i) 

2 - 

atmosphe r i c C0^ . 11 ) CC^ in carbonate minerals and iii) CO^ 
produced by oxidation of organic matter and carbon present in the 

rocks (Holland, 1978). As discussed in the earlier section, the 

second dissociation constant of H^CO^ has the value = 10 ^ ^ . 

It follows that at the pH value of river water which is much less 

2 - 

than 10.3, CO^ is relatively negligible compared to HCD^ . 
Therefore, among the dissolved carbonate species derived from the 
solubility of atmospheric CO^ gas. HCO^ is dominant in river 

water . 

The weakly acidic derived from atmospheric CO^ 

weathers carbonate and silicate minerals according to the 
following two representative reactions. 

a) Congruent dissolution of dolomite : 

CaMg ( CO j ^ + 2C0 2 + 2H 2 0 = Ca 2+ + Mg 2+ + 4HC0 ? _ 

Ca 2+ and Mg 2+ are balanced by HCO^ on an equivalent basis. 

b) In congruent decomposition of K-felspar : ^ 

2K-f elspar + 2CO £ + 11H 2 0 = Kaolinite + 2K + + ZHC0 ? + 4H 4 Si0 4 
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Therefore, the HCO^ in river water has two fractions viz. 

i) HCO^ derived from the weathering of primary carbonate 

minerals as well as calcareous cement. 

ii) HwOj derived from the weathering of silicate minerals. 

The total Ca and Mg in river water also have the same 
two sources, neglecting the weathering of evaporites. Hence the 
proportion of the carbonate and silicate fraction of total HC0~ 
can be calculated by the following steps suggested by Raymahashay 
< 1 986 > . 


1 . 

<HC0 5 \ O t 

SS 

(HC0 ? ) 

+ 

c 

CHCO,) . 

3 s i 

2. 

<HC0,> 

3 c 

* 

CCa) 

c 

+ 

(Mg) 

c 

3. 

<Ca, tot 

s 

CCa) 

c 

+ 

(Ca) . 

S 1 

4. 

<M9> tot 

a 

(Mg) 


(Mg) . 

S 1 

5. 

CHCO, > 

3 c 

a 

0.74<Ca) to t + 0 . 4(Mg) tot 

6. 

CHCO,) 

3 s i 

= 

(Hccy 

tot 

X 

o 

o 

V 

n 

where all concentrations 

are 

in mi 1 1 i -equ i val ents per litre and 

the 

subscr l pts 

s c’ 

, N S l f 

and 

'tot' stand for carbonate, silicate 


and total respectively. 

The results of this calculation for the Mahanadi river 

have been presented in Table 6.2. It is realized that the 

estimate of the proportion of HCO^ derived from carbonate 

weathering can be improved if factors for Ca and Mg specific to 

the Mahanadi basin are determined to replace the global average of 
0.74 and 0,40 adopted from Holland (1978). 

The ratio (HCOjW (HCCy g . serves as an index of 
carbonate versus silicate weathering in different parts of the 
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Table * 6.2 Bicarbonate derived from rock weathering in the 
Mahanadi basin C All concentrations are in 

milliequivalents per litre} 


Parameter C* tot Mg tot <HC0 3 > tot <HC0 3 > c <HC0 3 > si <HC0 3 > c /<HC 


Stn- 1 

0.65 

0.24 

0.84 

Stn-2 

0.91 

0.51 

1.21 

Stn- 5 

1 . 95 

0.87 

2.52 

Stn-4 

1 . 98 

1 .48 

5.64 

Stn- 5 

1 . 54 

1.04 

2.55 

Stn-6 

1 - 55 

0.55 

1.87 

Stn-7 

0.88 

0.58 

1.20 

Stn-8 

1 .17 

0.64 

1.81 

Stn-9 

0.57 

0.51 

0.91 

Stn- 1 0 

0.68 

0.55 

1.11 

Stn-11 

1 . 52 

0.66 

2.25 

Stn-12 

1.05 

0.52 

1.68 

Stn-2 5 

1 . 01 

0.41 

1.44 


0.56 

0.28 

2.00 

0.80 

0.41 

1.95 

1.49 

0.85 

1.79 

2.06 

1 . 58 

1.50 

1 . 56 

0.79 

1.97 

1.20 

0.67 

1 .79 

0.80 

0.40 

2.00 

1.12 

0.69 

1.62 

0.55 

0.56 

1.55 

0.64 

0.47 

1 .56 

1.24 

1 .01 

1 .25 

0.97 

0.71 

1.57 

0.91 

0.55 

1.72 
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basin. For example, the relatively high values of this ratio at 
stations 1, 2, 5 .5 and 6 probably reflect the dominance of 
carbonate weathering over silicate weathering. The geologic map 
(Fig. 6.2) and lithologic data (Table 6.3) support this 
conclusion. The high ratio at station 7 is at least partly due to 
the occurrence of carbonaceous shale in the drainage area. As 
discussed earlier, oxidation of carbonaceous matter in rocks is 
one of the sources of bicarbonate in river water. In contrast, 
the relatively low values at stations 9, 10, 11 and 12 may be due 
to the weathering of shale, sandstone, granite. khondalite. 
charnock i t e and similar silicate dominated rocks occurring in the 
drainage area. The low value at station 4 located within the 
limestone rich Chhattisgarh basin requires special attention. The 
main rock type in this area happens to be a reddish brown shale in 
contact with argillaceous dolomitic limestone (Das et al . , 1990). 
The X-ray pattern of this rock did not show any carbonate peaks. 
It appears that the weathering of the silicate rich shale 
surpasses the effect of carbonate weathering at this location. 
The (HCO, ) / ( HCO, ) . ratio at stations 8 and 13 has a value 
intermediate between the range of carbonate weathering and 
silicate weathering. These two stations are located on the main 
Mahanad i river channel. Therefore, it is highly probable that the 
river water composition at these stations reflect _ a combined 
influence of weathering of carbonate and silicate rocks which 
occur in the drainage basins of the tributaries. 
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Tabl© 6-3 Lithology around CWC sampling stations along the 
Mahanadi river 


Station Inferred rock types in the catchment area 


1- Limestone, granite and sandstone 

Z. Granite, limestone and sandstone 

3 . Dolomitic limestone, calcareous shale, sandstone and 

gran l te . 

4. Shale, dolomitic limestone and sandstone 

3. Dolomitic limestone, shale, sandstone and granite 

6. Granite, limestone and sandstone 

7. Carbonaceous shale, sandstone, granite and limestone 

8. Dolomitic limestone, shale, granite and sandstone 

9. Sandstone, shale and limestone 

10. Granite, sandstone and shale 

11. Granite and charnockite 

12. Granite, khondalite and charnockite 

1>. All the above rocks with alluvium. 


Ill 


6* 3* 3 Bsidric® of C 0«i + Mq) ^ridi HCQ 

3 

As discussed above, Ca, Mg and HCO ? in river water have 
two major sources namely, carbonate and silicate weathering. 

Therefore , 


(Ca + Mg> tot = (Ca + Mg> c + (Ca + Mg> si 

and <HC0 3 > tot = (HCO^ + (HCO^.. 

On the other hand, the stoichiometry of carbonate weathering 

reactions demands that (Ca + Mg) = (HCQ,) . However, the data 

c 3 c 

presented in Table 6.4 shows that in the Mahanadi river basin, 
(Ca + Mg) is not very different from (HCO, ) , t . Fig. 6.3 also 
shows that the water compositions fall very close to the equiline 
(Ca + Mg ) = (HCO^)^^. It is interesting to note that Stal lard 

and Edmond (1983) drew a similar conclusion for the Amazon basin. 


The near-equal i ty of (Ca + Mg) and (HCO,). , can be achieved in 

tot 3 tot 


the equation 


(Ca + Mg) + (Ca + Mg) . = (HCO,) + (HCO,) . 

c si 3c 3 si 

if the following two conditions are satisfied : 

l) (Ca + Mg) > (Ca + Mg) . and (HCO,) > (HCO, ) . 

C SI 3c 3 Si 

i . e . carbonate weathering dominates over silicate weathering. 

ii) (Ca + Mg) = (HCO, ) . 

SI /SI 

The estimation of carbonate weathering fraction of (Ca + Mg) and 
HCO^ from Holland’s global factor of (Ca> c = 0.74(Ca)^ ot and (Mg> c 
= 0 . 4 (Mg) tot indicates that the first condition is satisfied^ for 
each of the water samples of the Mahanadi river basin. In other 
words, carbonate weathering is relatively more important than 
silicate weathering throughout the Mahanadi basin. 
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Table 6.4 x Fractionation of Ca, Mg and HCC> 3 in water 


Parameter 

( Ca+Mg > 

tot 

<HC0, ) , 

3 tot 

(Ca+Mg) =<HC0 > (Ca+Mg> . 

c ? C SI 

<HCO t > 

3 s 

Stn-1 

0.87 

0.84 

0.56 

0.31 

0.28 

Stn-2 

1 .25 

1.21 

0.80 

0.43 

0.41 

Stn-3 

2.42 

2.32 

1.49 

0.93 

0.83 

Stn-4 

3.46 

3.64 

2.06 

1.40 

1.58 

Stn-5 

2.58 

2.35 

1.56 

1.02 

0.79 

Stn-6 

1 .88 

1 .87 

1.20 

0.68 

0.67 

Stn-7 

1 .26 

1 .20 

0.80 

0.46 

0.40 

Stn-8 

1.82 

1.81 

1.13 

0.69 

0.68 

Stn-9 

0.88 

0.91 

0.55 

0.33 

0.36 

Stn-1 0 

1.01 

1.11 

0.64 

0.37 

0.47 

Stn-1 1 

1 .98 

2.25 

1.24 

0.74 

1.01 

Stn-12 

1 .55 

1.68 

0.97 

0 . 58 

0.71 

Stn-1 5 

1 . 42 

1.44 

0.91 

0 .51 

0.53 


(All concentrations are in mi 1 1 i equ i val ents per litre) 



(Ca + Mg) tot meq/ 
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(HC0 3 )tot meq/l 


Fig. 6.2 Balance of CCa+Mgl total and C HCO^ total 
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There are two local factors which might explain the 
small deviations from the equiline in Fig. 6.3. For example, the 
weathering of evaporite beds has been neglected in these 
calculations. However, it has been realized that the small excess 
of <Ca + over (HCO^>^_ o ^ at stations 3 and 5 can be produced 
by this extra source of Ca and Mg. In fact gypsum occurs as 
cementing material in surface exposures of shale and mudstone in 
the Seonath river section. It has been observed in bore hole 
sections within calcareous and non -cal careous shales of this area 
(Das et al., 1990). Similarly, the small deficiency in (Ca + 

at stations 4, 10, 11 and 12 can be produced by 

2 + 

ion-exchange between Na-clays and Ca in water (Cerling et al . , 
1989). This mechan i sm is feasible in the silicate dominated 
drainage areas corresponding to these stations. 


6.3.4 Ternary Diagrams 

The relative role of carbonate, silicate and evaporite 

weathering in influencing river water composition can be evaluated 

through a ternary diagram (Fig, 6.3) similar to that developed by 

Stal lard and Edmonds (1983). The relevant percentages of HCQ V 

S l Q_ and SO, are given in Table 6.3. In the diagram, the Mahanadi 
2 4 

river waters plot near the bicarbonate vertex. This indicates 
that the weathering of carbonate rocks has the maximum influence 
on river water composition. 

The charge balance of (Ca 2+ + Mg 2+ > tot by HC0 ? or SO 2 
can be further verified by another ternary plot of these 
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Table 6.5 Relative percentages of HCO , SiO and 

3 2 

an equivalent basis 


SO A on 

4 


Samp 1 i ng 
stations 

hco 3 

Si0 2 

S0 4 

1 . 

63.57 

18.19 

16.24 

2. 

75.95 

14.90 

9.15 

3. 

77.48 

7.84 

14.68 

4. 

70.63 

6.18 

23.19 

5. 

71 .29 

8.28 

20.43 

6 . 

80.52 

10.19 

9.29 

7. 

64.40 

13.16 

22.79 

8 . 

72.95 

11.14 

15.91 

9. 

59.91 

18.96 

21.13 

10 . 

67.23 

16.05 

16.72 

11 . 

83.92 

9.08 

7.00 

12. 

76.34 

13.81 

9.85 

13. 

76.33 

12.82 

10.85 
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parameters (Fig. 6.4). The relevant percentages are given in 

Table 6.6. It is clearly seen that water compositions from the 

2 + ? + — 

Mahanadi river fall along the line (Ca + Mg ) = HCD^ which once 
again confirms that weathering of carbonates is the dominant 

process. 

6.4 Carbonate Mineral Equilibria 

The broad geochemical trends discussed in the earlier 
sections indicate that weathering of dolomitic limestone has a 
pronounced influence on water chemistry in the upper part of the 
Mahanadi river basin. Therefore, it would be important to 
evaluate the state of saturation of river water in terms of the 
solubility products of the two minerals calcite and dolomite. The 
relevant reactions are 


i) CaCO. 


K 


a Ca 2 + 


i i > CaMg< CO ^ 
K , 3 


2 + 2 - 
Ca 4 + CO^ 

a. Z- = lCf 8 ' 4 at 25°C. 
CO, 

? 

* Ca 2 + + Mg 2 + + 2C0 2 


a Ca 2 + 


a Mg 2 + 


(a 2- ) 2 = 10“ 17 -° at 25°C. 

va co: 

> 


At the first instance, the ionic activity product UAP> 
in river wate r was calculated using pH, ( ^ tot* Mg tot 


(a 


HCO ? tot 


1 isted in Table 6.7. The activity of CO, 


ion was 


, , , , - /_ _> from the relationship 

calculated from <a nco^ ; tot 
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Table 6*6 

Relative percentages of CCa+MgD, 

equivalent basis 

HCO and SO on 

o 4 

Sampl i ng 

(Ca+Mg) 

HC0 ? 

S0 4 

1 . 

45.47 

45.71 

10.82 

2. 

47.51 

46.85 

5.64 

5. 

46.67 

44.84 

8.49 

4. 

41 .71 

44.17 

14.52 

5. 

46.08 

41.92 

12.00 

6. 

47.41 

47.15 

5 . 44 

7. 

45.70 

41.52 

14.78 

8. 

45.24 

44.96 

9.80 

9. 

41.99 

42.88 

15.15 

10 . 

42.50 

46.21 

11.49 

11 . 

44.75 

51.00 

4.25 

12. 

45.14 

48.59 

6.27 

n. 

46.59 

46.94 

6.67 
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a CQ^" ‘ a H + 

k 2 = = 10" 10 - 3 at 25°C . 

^ a HC0~ } tot 

The calculated values of IAP are also listed in Table 6.7. It is 
clear that the river water at all stations is undersaturated with 
respect to calcite and dolomite with the exception of station 4. 
This should account for the aggressive weathering of dolomitic 
limestone throughout the river basin. 

An attempt was made to refine this calculation by 
considering the fraction of Ca and Mg derived from carbonate 
weathering instead of total Ca and Mg. As before on an equivalent 
bas i s 


< Ca > * 

c 

(Mg) = 
c 

Therefore. (HCO,) = 

> c 


0.74<Ca>. . 

tot 

0.4(Mg) tot 
< Ca) + (Mg) 


and 


a C0*~ 


K Z ' ^HCO^c 

V 


The revised calculations are presented in Table 6.8. 
This revision confirms the earlier conclusion that except at 
station 4, the river water is undersaturated with respect to 
calcite and dolomite throughout the basin. This trend is 
presented graphically in Figs. 6.5(a) and 6.5(b) where the 
logarithmic values of IAP for calcite and dolomite are plotted 
against TDS . It is obvious that at high values of dissolved ion 
concentration, the water approaches saturation with calcite and 


I 



carbonate derl vod Ca , Mg and HCO in water 




BO, [-H. +2 *3*1 Bo, ’°\ + 2 6 ^>B 0 , ^< + 2 ^>Bo, 
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dolomite. For example, TDS is maximum at station 4 and this may 
be the reason why the river water is practically saturated with 
respect to calcite and dolomite at this location. 

Similar conclusions can be drawn by evaluating a 

Saturation Index (SI> defined as 



It is obvious that a value of SI = 0 indicates saturation. The SI 
values have been presented in Table 6.9 and the same plotted in 
Fig. 6.6. The general undersaturation of the river water 

simultaneously with respect to calcite and dolomite is clear from 

this f i gure . 

The calcium and magnesium concentration in river water 
can also be used to evaluate the possibility of dolomit izat ion of 
calcite by limestone-water interaction. Calcite can be converted 

to dolomite by the reaction 

2CaC0 ? + Mg 2+ = CaMg (C0 3 > 2 + Ca 2 

The equilibrium constant of this reaction is expressed as 

^cal-dol 
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Table 6. Q Estimation of Saturation Index CSID of calcite and 

dolomite 


Stat i on 


SI 


calcite 


SI 


dolomi 


te 


1 . 

- 1.58 

- 3.64 

2 . 

- 1.29 

- 3.12 

3 . 

- 0.28 

- 0.89 

4 . 

0.00 

- 0.19 

5 . 

- 0.25 

- 0.74 

6 . 

- 0.34 

- 1.13 

7 . 

- 1.18 

- 2.79 

8 . 

- 0.51 

- 1.35 

9 . 

- 1.49 

- 3.30 

10 . 

- 1.28 

- 2.90 

11 . 

- 0.28 

- 0.92 

12 . 

- 0.68 

- 1.73 

13 . 

- 0.93 

- 2.31 


dolomite 
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Dolomite saturation 


GO 



• 4.0 -p“T T ""I r'T-T'-T 

- 1.7 - 1.2 


i i i i ii i i i 

0.7 


"i n i i f r r 

- 0.2 0-0 


SI 


caicite 


Fig. 6.6 Saturation indices for caicite and dolomite with 
respect to river water samples 


Caicite saturation 
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This is numerically equal to 

K* , ln -8.4 N 2 

- <10 > ,-0.20 

£r- - r r? ' /) — - 10 

d <10 1/ ’°> 

In other words calcite can be converted to dolomite only if the 
value of a Q a ^ + / / /a ^jg2 + is less than 10^‘ 2 ^. The value of this 
ratio at various stations along the Mahanadi river have been 
listed in Table 6.10. It is clear that the river water is not 
capable of do 1 omi 1 1 zat l on except at station 4 and 5. 

6.5 Thermodynamics of Silicate-Water Equilibria 

Thermodynamic data allow us to verify whether a given 
water environment is suitable for stability of silicate minerals 
found in weathered rocks, soils and sediments. 

To start with, we need a set of weathering reactions 
involving primary minerals, secondary minerals and dissolved ions. 
It is assumed that these reactions approach equilibrium in the 
weathering environment. In that case, the equilibrium constant 
(K) of these reactions can be estimated as discussed below. 

For a reversible reaction such as 
bB + cC ^ dD + eE 

a Q d . a^e 
K " a B b . a c cT 

where the symbol 'a* stands for activity or effective 
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Tafoi^ 6*10 Calculation f or dol onti ti zat i on of* calcito 


103 [ <a Ca 2 +> tot /<a M g 2 +, t at] 


1 . 

0.42 

2 . 

0.46 

5 . 

0.25 

4 . 

0.05 

5 . 

0.17 

6 . 

0.58 

7 . 

0.28 

8 . 

0.26 

9 . 

0.26 

10 . 

0.51 

11 , 

0.50 

12 . 

0.29 

15 . 

0.59 
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concentration of the species taking part in the reaction. For 
pure minerals and liquid H 2 0, activity is assumed to be equal to 

unity. 

This constant is related to the change in Gibbs free 
energy during reaction (AG^) under standard conditions (298 . 1 6°K 
and 1 atmosphere pressure) through the equation 



AG = 
r 

-RT In 

K 


or 

AG° = 
r 

-2.505 

RT log K 

(1 ) 


where R is universal gas constant equal to 1.98 x 10 

-3 0 

kca 1 / deg . mo 1 e or 8.28 x 10 kJ/deg.mole. AG^ for a given 
weathering reaction can be calculated from the relationship 

AG° = GAG° - I AG° 

r 'product 'reactant 

The symbol AG° stands for standard state Gibbs free energy of 
formation of various species participating in the reaction. 

At 2 5°C , the above equation (1) reduces to 

AG° * -1.564 log K 
r 

in kcal 

or AG° * - 5.707 log K 

r 


i n k3 
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6.5.1 Construction of Stability Diagrams 

Sufficient information are now available for the 
construction of stability diagrams in the four systems : 


a > 

1 

o 

Csl 

b> 

Na 0- 

c ) 

CaO - 

d> 

MgQ - 


M 2°5 - 

SiQ 2 

1 

X 

ho 

o 

A1 2°5 - 

Si0 2 

- H 2 0 

A1 2°> - 

Si0 2 

- H 2 0 

> 
t • 

0 

1 

Si0 2 

- H z 0 


As illustrated below, equilibrium constants of relevant weathering 
reactions are simple functions of cation to hydrogen ion activity 
ratio and silica activity. By convention, these weathering 
reactions are written with the 'conservation of A1 principle’ i.e 
all A 1 is locked up in the weathering products and the solution 
does not contain dissolved A1 . Thus stability diagrams can be 
constructed with logarithmic axes ashowing a^+^Z a^+ , a Na + / / / a H + ’ 

a Ca 2+ / <a H +>2 ’ a Mg 2+ / <a H +> and "V^SiO^ 

Three constant silica lines are common to the four 
diagrams < F i g . 6.7, 6.8, 6.9 and 6.10). These are 

1) Solubility of quartz . 6 ppm SiO^ at 25°C 

Quartz + 2H,,0 = H^SiO^ 

log Kj * log a H SiO 


- 4.0 
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2 ) Sol ub i 1 1 1 y of amorphous s i 1 i ca. , 1 


20 ppm Si0 2 at 25 C 


Si0 2 + 2H 2 0 = H 4 SiD 4 

log K. = log a = -Z 7 

2 H.SiO, t - m/ 

4 4 


5 > Gibbs 1 1 e-Kaol in i te boundary 

A1 2 Sl 2°5 <OH> 4 + ?H 2° = A1 2 (OH) ^ + 2H 4 Si0 4 

J/2 log K, = log a 

4 4 


The value of silica activity can be calculated from AG° 

r 

of this reaction. However, it has been observed (Drever, 1988) 

that large uncertainty prevails in the reported values of free 

energy of silicate minerals involved in weathering reactions. 

Therefore, those thermodynamic values have been selected which 

match with equilibrium constants independently estimated from 

experimental works and field observations. For example, accepting 

quartz solubility of 6 ppm and AG° of quartz and H 2 0 as -856.29 

and 257.14 kJ/mole (Robie et al . . 1978) respectively, AG° of 

H.SiO. can be calculated as follows : 

4 4 

Quartz + 2H 2 0 = H 4 Si0 4 

log K = - log a R SiQ = - 4.0 

4 4 

AG° * 4 X 5.707 = 22.85 kJoules 

r 

AG° H 4 Si0 4 * - 1507.74 kJ/mole. 
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This value differs from - 1316.6 kj/mole used by Drever 

(1988) to construct the stability diagram. On the other hand, it 

matches closely with majority of values listed by Woods and 

Garre i s ( 1 987) including -1308.0 reported by Robie et al . (1978). 

Using this value of AG° H 4 Si0 4 and-2309.78 kJ/mole for 

gibbsite -3799.36 kj/mole for kaolinite (Robie et al . . 1978), AG° 

r 

for the g i bbs i te-kao 1 i n i te reaction is 59.80. 

j log = — (4|77iy7> = -5.24 

_ t .5 ?4 

This value of a u e . n = 10 ' for g i bbs i te-kao 1 i n i te 

H 4 5l0 4 

-4 4 

equilibrium obviously differs from 10 ' used by Drever (1988) 

but is close to 10 ^ derived by Krauskopf (1979) from a different 

set of f ree-energy values. 

K 2 0_A1 2°^" Si °2~ H 2° SYSTEM 

The boundaries of the stability diagrams for the 
K-system have been drawn by calculating equilibrium constants for 

the following reactions. 

4) Muscov i te-Kaol i n i te boundary 


2KAl 3 SijO 10 (OH> 2 + 2H + + 3H £ 0 


3Al,Si,0e<0H) . + 2K + 
2 2? 4 


From the data in Table 6.11 

AG° * -50.3 kJoules 

r 



Table 6.11 Free energy values CkJy'moleD for 
solutes at 25 C temperature and 1 

var i ous mi ner al s and 

atmosphere 

Gibbsite CAl.(OH)^] 

Z 6 

-2309.78 

Kaol mite [Al S i ,0. (OH) . ] 

2 2 5 4 

-3799.36 

Muscovite [KA1 ^Si ^0 1Q (0H) 2 3 

-5600.67 

K- f e 1 spar [KAlSi.O.] 

> 8 

-3742.33 

Na- smect i te [Na Q ?? A1 2 ^Si^ 67 ° 10 <OH) 2 ] 

-5368.39 <;i> 

Ca-smect i te CCa Q 16/ A1 2 ^^Si ? 67 ° 10 <OH) 2 :i 

< 2 ) 

-5376.96 

Chlorite tM 9 5 Ai 2 S ! 5 ° 10 <G H > 8 3 

-8207. 75 0> 

Quartz CS.0.,3 

- 856.29 

Na + 

- 261.90 

K + 

- 282.49 

Ca 2 + 

- 553.54 

.. 2 + 

Mg 

- 454.80 

H , S i 0 . 

4 4 

-1307.74 <4) 

H n 

- 237.14 

2 



<1> Calculated by extrapolation of Hemley's data (Hess, 1966) 

(2) Calculated from log K = -18.4 for Ca smect i te-kao 1 i n i te 

equilibrium (Norton, 1974) 

(3) Helgeson et al , (1978) 

Calculated from 6 ppm SiO^ for quartz solubility 
All other values from Robie et al (1978). 


( 4 ) 
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K 2 0 - AI2O3 - S102 " H 2O System 


FIG 6-7 Stability diagram: 
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log « 4 = 2 log <a K +/a H +) =8.81 

log <a K +/a H +> = 4.40 

5 )Muscovite-Gi bbs i te boundary 


2KA1 3 Si 3 O 10 <OH) 2 + 2H + + 18H 2 0 = 3A1 2 <0H> 6 + 2K + + 6H 4 Si0 4 


A line of slope -3:1 passes through the point of 

intersection of line (3) and (4). 


6) K- f e 1 spar-Kao 1 i n i te boundary 


2KA1 S i ,0 D + 2H + + 9H,0 = Al o Si„0 t .<0H). + 2K + + 4H.SiO / 
38 2 2254 44 

AG° = 23.62 kJoules 

r 

log = -4.14 

1/2 log K 6 = log (a K +/a H +> + 2 log a R SiQ = -2.07 


7> K- f © 1 spar-Muscov i te boundary 


3KA1 S i ^D g + 2H + + 12H 2 0 = KA1 i ^ Q < OH) £ + 6H 4 Si0 4 + 2K + 

A line of slope -3:1 passes through the point of 
intersection of (4) and (6). 


Na 2 0-Al 2 0 3 -Si0 2 -H 2 0 SYSTEM 

Stab i 1 i ty d i agram of the Na — s y stem con ta ins a new line 

representing equilibrium between kaolinite and an ideal 

1 HI 



Log ( a Na V a ^ ) 
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FIG 6-8 Stability diagram: Na20 -AI2O3 - Si O2 - H2O System 
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Na-montmor i 1 loni te 


with 


structural 


formula 


Na 0 . 33 Ai 2 . 33^ * 3 . 67^1 0^ HJ 2 ' the other hand, there is a great 
deal of confusion about the AG value of this phase as well as the 
exact composition of montmor i 1 Ion i te formed in weathering 
environment. Some of these aspects have been discussed by Drever 
(1988). Therefore, the equilibrium constant for kaolinite-Na 
montmor i 1 1 on i te reaction was estimated in the following way. 


8) Na montmor i 1 loni te-Kaol ini te boundary 

6Na Q ^Al 2 3 ^Si 3 67 0 1() (0H> 2 + 2H + + 23H 2 0 

* 7A1 -S i o 0, (OH) , + 8H. S iO + 2Na + 

L L ? 4 4 4 

log K 8 = 21og <a Na +/ a H +> + 81og a 

4 4 

Hess (1966) extrapolated. Hemley’s high temperature 
experimental data to 25°C and estimated log (a Na +/a H +> = 5.9 and 

log a „ n * -3.5 for the above equilibrium. 

H 4 SlD 4 

.*. log K 0 * -14.6 

The equation for Na montmori 1 lonite-kaol ini te boundary 

i s 

log <» Na + / a H +> + 4109 a H 4 Si0 4 = 1/2 losK 8 = ~ 7 * 5 

AG° = 83.32 kJoules. Knowing the AG° of other 

r 

specie. (Table 6.11). AG° of Na-montmor i Uonite corresponding to 
the above equilibrium constant can be calculated as -5568.59 


kJ/mol e . 
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Table 

1 ) 

2 ) 

3> 

4) 

5> 

6) 

7> 


Values of AG^ of Na-montmorillonite as reported in 
recent literature 


Raymahashay <1968) 

He 1 geson (1969) 

Nnagu (1975) 

Woods & Garrels (1987) 
Drever (1988) 
Karathanas l s (1989) 

This report a 


-1277.4 kcal/mole 

— 5344.64 kJ/mole 

-1278.9 kcal/mole 

-5350.92 kJ/mole 

-1277.76 kcal/mole 

-5350.92 kO/mole 

-5346.1 kJ/mole 

-5382.0 kJ/mole 

-5080.5 to -5262.0 kJ/mole 
(soil smectite) 

-5368.39 kJ/mole 


3 Calculated by extrapolation of Hemley’s data (Hess. 1966) 
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This value has been compared with those reported in 
recent literature (Table 6.12). In the resulting stability 
diagram (Fig. 6.8), the kaolinite field is slightly wider compared 
with the diagram recommended by Drever (1988) 

Ca0-AI 2 0 3 -Si0 2 -H 2 0 SYSTEM 

Stability diagram for the calcium system contains the 
boundary between kaolinite and an ideal Ca- montmor ill on i te 

^ a 0 . 1 67 ^* 2 33^*3 67 ^ 10 ^*~^ 2 ^ * n manner similar to the sodium 

system. The equilibrium consant for Ca montmor i 1 1 on i te-kao 1 i n i te 

reaction was calculated in the following way. 

9 ) Ca montmor i 1 1 on i te-Kaol i ni te boundary 

6Ca 0.167 A1 2.5> S h.67°10 (DH> 2 + 2H+ + ” H 2° 

= 7Al 2 Si 2 CyOH> 4 + 8H 4 Si0 4 + Ca 2+ 

log K g * log [a Ca 2+/ <a H +> 2 ] + 81 ° 9 a H 4 Si0 4 

From his study of river water and sediments in the 
Rio Tanama system, Puerto Rico, Norton (1974), calculated 
the value for log K as - 18.4 ± 0.8. Therefore, the 

Ca montmor i 1 loni te-kaol inite boundary can have the equation 

log [a Ca 2V<V >Z ] + 8109 V i0 4 = ' 18 '‘ 
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Robie etal (1978) 

— Drever (1988 ) 



FIG 6.9 Stabili ty diagram : CaO - A^C^-SiC^ - H 2 O System 
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Table 6.13 Values of AG° Ca-montmorillonite as reported in 
recent literature 

1. Heigeson <1969) s -1279.24 kcal/mole 

= -5352.34 kJ/mole 

2) Nr i agu <1975) : -1280.9 kcal/mole 

= -5359.28 kJ/mole 

*1 

3> Raymahashay (1984) : -1278.84 kcal/mole 

= -5350.67 kJ/mole 

4) Woods & Garre 1 s (1987) : -5349.1 to -5354.9 kJ/mole 

5) Drever <1988) : -5388.0 kJ/mole 

6) Karathanasis (1989) : -5374.0 kJ/mole 

(soil smectite) 

7) This report : -5376.96 kJ/mole 

*1 Calculated from log K * -15.7 for Ca montmor i 1 Ion i te-kaol ini te 

equilibrium (Tardy. 1971). 

*2 Calculated from log K = -18.4 for Ca-montmor i 1 1 on i te-kao 1 i ni te 

equilibrium (Norton, 1974). 
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Using the above log K, 

AG° = 105.0 kJoules. 

Knowing the AG f of other species (Table 6.11), AG° of 
Ca-montmori 1 loni te can be calculated as -5576.96 kJ/mole. This 
value has been compared with other sources in Table 6.15. As is 
seen from the resulting stability diagram (Fig. 6.11), the 
kaolinite field has slightly narrower width compared with that 
recommended by Drever (1988). 

MgO-Al 2 0 3 -Si 0 2 -H 2 0 SYSTEM 

Stability diagram for the magnesium system shows two new 
boundaries i) between kaolinite and chlorite, Mg^Al ^Si ^0 ^ ^ (OH) g 
and ii) between gibbsite and chlorite. The diagram has been 

constructed in the following way. 

1 0 ) Chlor i te-Kaol ini te boundary 

Mg 5 Al 2 Si 3 O 10 (OH) 8 + 1 0H + = A1 2 S i ^ (OH) 4 + 5Mg 2+ + H 4 Si0 4 + 5H 2 0 

From the AG° values listed in the Table 6.11, 

AG° = -559.05 kJoules 

r 

log K 1Q = 62.91 

5 log K = log [a Mg 2 + / <a H +> 2 ] + 109 a H 4 Si0 4 = 


or 


62.91 
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FIG 6.10 Stability diagram.-MgO - Al 2 0 3 -Si0 2 - H 2 0 System 
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[1° 9 *h siO. 

4 4 


This line intersects the g i bbs i te-kao 1 i n i te boundary 
-5.24] at log [> Mg 2+/(a H + > 2 ] = 13.63. 


11) Ch 1 or l te-G i bbsj^te boundary 


Mg 5 Al 2 Si ? O 10 <OH> 8 + 10H T = A1 (OH) + 3Mg 2 + + 3H SiO 


A line of slope -3:5 through the point of intersection 

of g i bbs 1 1 e -kao 1 l n i t e and ch 1 or i te-kao 1 i n i te boundaries. The 

stability diagram (Fig. 6.12) also shows the relative position of 

chlorite field based on gibbsi te-kaol ini te boundary at log a H S . Q 

4 4 

* 4.4 (Drever . 1988). 

It is clear from the stability diagrams that positions 
of the boundaries in the K, Na. Ca and Mg system are slightly 
different from the diagrams constructed by Drever (1988). It is 
to be recalled that Drever selected AG° values to match the 
boundaries with natural water composition which were expected to 
be in equilibrium with the appropriate minerals. On the other 
hand, the boundaries in the stability diagrams are based on 
thermodynamic data mostly from Robie et al . (1978). Another 

reason for the difference between the two sets of diagrams is the 
AG° value of H 4 Si0 4 calculated in this work (-1307.74 kJ/mole> 
from 6 ppm SiO^ for quartz solubility which is higher than the 
value -1316.6 kJ/mole used by Drever (1988) to construct his 
stability diagrams. However, a value of -1308.1 kJ/mole has been 
listed in the appendix of Drever (1988) after adjustment for 
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consistency with other values, 
used in this thesis. 


This is almost equal to the value 


6. 5. 2 Silicate Hi ner al “Water Equilibria 

Once the stability diagrams were constructed as 
discussed above, the river water chemistry was evaluated by 
calculating the relevant ion activities from concentrations in 
mole/1 itre, ionic strength and activity co-efficients. The sodium 
concentration was corrected for cyclic salt contribution by 
assuming that 0.39 ppm chloride in rainwater collected at Raipur, 
came as NaCl. Therefore, an equivalent amount of sodium (0.017 
mi 1 1 l mo 1 e / 1 t t re ) was subtracted from the sodium value in river 
water. The appropriate calculations for plotting river water 
analyses on the stability diagrams have been summarized in Tables 
6.14, 6.13, 6.16 and 6.17. The position of water compositions 
show the following features. 

i> Kaolinite is a stable mineral in the river water 
environment with reference to potassium, sodium, calcium and 
magnesium systems. This is supported by the observation that this 
mineral is a common component of suspended sediments and bank 
deposits in the samples collected from the upper part of the basin 
during this investigation. Naidu et al . (1985) have reported that 
kaolinite is a constituent of the bed sediments in the downstream 
part of river. Similarly. Chakrapani and Subramanian (1990a) 
also found kaolinite in both bed and suspended sediments at 
various locations in the Mahanadi river basin. 
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Table 6.14 Water composition data used in the stability diagram 
for K a O-Al 2 0 3 -Si0 2 -H 2 0 system 


Stn. 

pH 

K + 

m . mo 1 e / 1 it 

a K + . 10 3 

log(a K +/a H +) 

log a H 

1) 

7.23 

0.052 

0.050 

2.929 

-3.633 

2) 

7.21 

0.050 

0.047 

2.882 

-3.623 

3) 

7.73 

0.128 

0.118 

3.822 

-3.629 

4) 

7.76 

0.106 

0.097 

3.745 

-3.496 

5) 

7.77 

0.100 

0.092 

3.734 

-3.564 

6) 

7.57 

0.048 

0.045 

3.223 

-3.625 

7) 

7.35 

0.078 

0.074 

3.219 

-3.609 

8) 

7.73 

0.080 

0.075 

3.605 

-3.559 

9) 

7.33 

0.057 

0.054 

3.062 

-3.542 

10) 

7.39 

0.038 

0.036 

2.946 

-3.577 

11) 

7.83 

0.045 

0.042 

3.453 

-3.613 

12) 

7 . 64 

0.056 

0.053 

3.364 

-3.519 

13) 

7.47 

0.044 

0.042 

3-113 

-3.618 
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Table 6* 1 5 


Water composition data used in the stability diagram 
for Na a 0-Al 2 0 3 -Si0 2 -H 2 0 system 


Stn. 

pH 

Na + 

m . mo 1 e / 1 it 

Cyl 

Na + 

m. mole/1 it 

ic salt 

W- 10 ’ 

corrected* 

lo g(a N a+/a H +) 

109 a H,S i 0. 

4 4 

1. 

7.23 

0.210 

0.193 

0.184 

3.494 

-3.633 

2. 

7.21 

0 . 307 

0.290 

0.269 

3.639 

-3.623 

3. 

7.75 

0.852 

0.835 

0.772 

4.638 

-3.629 

4. 

7.76 

1 .896 

1.879 

1 .713 

4.994 

-3.496 

5. 

7.77 

0.989 

0.972 

0.899 

4.724 

-3.564 

6) 

7.57 

0.592 

0.575 

0.538 

4.301 

-3.625 

7) 

7.35 

0.486 

0.469 

0.442 

3.995 

-3.609 

8) 

7.73 

0.674 

0.657 

0.615 

4.519 

-3.559 

9) 

7.33 

0 . 248 

0.231 

0.220 

3.672 

-3 . 542 

10) 

7.39 

0 .424 

0.407 

0.386 

3.976 

-3.577 

11) 

7.83 

0.970 

0.953 

0.890 

4.779 

-3.613 

12) 

7.64 

0.643 

0.626 

0.589 

4.410 

-3.519 

13) 

7.47 

0.342 

0.323 

0.307 

3.957 

-3.613 


* 0.39 ppm or 0.017 mi 1 1 i -mol e/ 1 i tre of chloride equ i val ent-Na + 
in rain water has been subtracted from the total sodium 
concentration in river water. 
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Table 6.16 

Water composition i 

for Ca0-Al 2 0 3 ~Si0 2 

data 

- H 2° 

used in 

system 

the 

stability 

diagram 

Stn . 

pH <Ca 2+ ) tot 

m.mola/1 it 

< a r 2 + ) 

Ca tot 

xl0 ? 

log 


tot 

r(a 2+) . 

loal Ca Sl 

] 109 a H 4 SiO 

L(a H+ ) 2 -i 

9 2 

L <a H +) 

1) 

7.23 

0.317 

0.262 


10.878 


10.293 

-3.633 

2 > 

7.21 

0.458 

0.370 


10.988 


10.403 

-3.623 

3) 

7.73 

0.776 

0.567 


12.254 


11.669 

-3.629 

4) 

7.76 

0.989 

0.683 


12.354 


11.769 

-3.496 

5) 

7.77 

0.771 

0.564 


12.291 


11 .706 

-3.564 

6) 

7. 57 

0.665 

0.512 


11.849 


11.164 

-3.625 

7) 

7.35 

0.441 

0.349 


11.243 


10.658 

-3.609 

8 > 

7.73 

0 . 587 

0.450 


12.113 


11.528 

-3.559 

9) 

7.33 

0.287 

0.235 


11.031 


10.446 

-3.542 

10) 

7.39 

0.341 

0.276 


11.220 


10.634 

-3.577 

11) 

7.83 

0.660 

0.501 


12.360 


11.775 

-3.613 

12) 

7.64 

0.516 

0.404 


11.886 


11.301 

-3.519 

13) 

7.47 

0.504 

0.399 


11.541 


10.956 

-3.618 
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Table 6.17 Water composition data used in the stability diagram 
for MgO-AlgO^-SiOg-HLgO system 


Stn. 

pH 

2 + 

<M9 ’tot 

m . mo 1 e / 1 it 

<a.. 2 + ) , , 

Mg tot 

xl0 ? 

log fW + 1 

109 a H SiO 
4 

La H +) 2 Jtot 

1) 

7.23 

0.120 

0.099 

10.456 

-3.633 

2) 

7.21 

0.157 

0.127 

10.524 

-3.623 

5) 

7.75 

0.433 

0.317 

12.001 

-3.629 

4) 

7.76 

0.742 

0.512 

12.229 

-3.496 

5) 

7.77 

0.521 

0.381 

12.121 

-3.564 

6) 

7.57 

0.276 

0.212 

11 .466 

-3.625 

7) 

7.35 

0.189 

0.150 

10.876 

-3.609 

8) 

7.73 

0.322 

0.247 

11.853 

-3.559 

9> 

7.33 

0.157 

0.129 

10.771 

-3.542 

10> 

7.39 

0.167 

0.135 

10.910 

-3.577 

11) 

7.83 

0.329 

0.250 

12.058 

-3.613 

12) 

7.64 

0.262 

0.205 

11.592 

-3.519 

13) 

7.47 

0.205 

0.162 

11.150 

-3.618 
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ii> In the stability diagram for Ca-system <Fig. 6.9), 
the water compositions plot near the kaolinite-Ca montmori 1 Ionite 
boundary. This indicates that the kaolinite-Ca montmori 1 Ionite 
equilibrium controls river water composition. A correction for Ca 
concentration was introduced by considering the fraction derived 
from weathering of silicate minerals instead of the total Ca 
value. As discussed earlier 


<Ca> 


(Ca). - (Ca) 

tot c 

1 - 0.74<CaK . 


= 0.26 (Ca). . 

tot 

The corrected values of a Ca 2+ / (a H +> | have been listed in the 

L s i 

penultimate column of Table 6.16. The shift of position of the 
stability fields is negligible and the earlier conclusion that the 
kaolinite-Ca montmor i 1 1 on i te equilibrium might control river water 
composition, remains valid. This idea appears to be feasible 
because suspended sediments and soil samples from Andhiyarkore 
(station 4) as well as soil samples from Simga (station 3) and 
Jondhra (station 5) showed a mixture of kaolinite and glycol 
expansive clay minerals during X-ray diffraction studies. 
However, it was not possible to confirm the exact composition of 
montmor i 1 1 on i te . Previous workers have reported montmor i 1 1 on i te 


in the suspended and bad loads of the Mahanadi river particularly 
in the lower reaches (Naidu et al . . 1985; Chakrapani and 

Subramaman. 1990a). There is a strong possibility therefore, 
that a calcium rich montmor i 1 Ion i te forms by weathering of 
limestone and calcareous shale in the drainage basin. 



6.6 Aluminiuro Spec! at ion 


The stability diagrams developed in the earlier section 
were based on the assumption that aluminium is immobile during 
rock weathering and does not appear in the solution phase. On the 
other hand, the range of A1 concentration listed in Table 5.1 
indicates that the maximum concentration reaches nearly 1 ppm. 
Therefore, it is important to verify if the stability of an 
aluminium bearing mineral has any control on dissolved A1 at this 
concentration level . Since kaoiinite is al ways present as an 
Al-siiicate phase in suspended sediments and soils within the 
Mahanad l basin, its solubility in terms of dissolved Al species 
has been evaluated in this section. 

It is well established that the major dissolved Al 

3+ 2+ + 0 

species in natural waters are Al , Al (OH) , A 1 < OH ) 2 , A 1 < OH ) ^ and 
Al (OH) . The total dissolved aluminium in moles per litre which 
is the analytical value, can be expressed as 


m Al T = m Al 3 * * m Al (OH) 2 * + m Al<OH); + m AKOH>; + m Al<QH>; (A> 

The activities of these dissolved species are related to one 
another through the following equilibrium constants. 

Al 5+ + H -0 = A 1 ( 0 H) 2+ + H + ( 1 ) 

£ 

a Al<OH> 2+ * a H + 
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Therefore 


S irru lar 1 y 


r Z • "’aHOH) 2 ' 1 ' ■ a H* 
^5 ' m Al ,+ 


K i • r -i • m A1 3+ 

m Al(OH) 2+ = 

Y 2 ’ m Al 5+ 


Al ,+ + 2 H 2 0 


A 1 < 0 H>* + 2 H + 


( 2 ) 


K * 
2 


a Al(OH > 2 • <a H +> 
a Al 3 + 


Y \ • m Al(OH>: • <a H +) ' 


- m A1 ? + 


m A 1 < OH ) 2 


K 2 • h • m Al 3+ 
r l * (a H +>2 


Al^ + + = A 1 ( 0 H)° + 5 H H 


( 3 ) 


K_ = 
> 


a Al(OH)° • <a H +> 

i~Tf 

AI 


• m A l < OH ) ° ■ <a H +> 
' m Al ,+ 
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m Al (0H>° 


K 3 * • m A ,3+ 


A1 


‘ <V> 


C4> 


A1 


5+ 


+ 4H £ 0 = A1 (OH) ~ + 4H + 


K, 


AKQH> 4 . <a H +)‘ 
a Al ?+ 


r A ' m Al(QH)' * <a H +> 
4 

y 3 • m Al 5t 


m 


A 1 < OH > 


K 4 ' Y J_ • m Al ? + 
r 4 * <a H +> 


Substituting these values in equation (4), 


'V ) V m Al ,+ V'V’aI 5 * K v>V m A1 3+ 

m Al = m Al 5+ + 2 3 

T y , . a H + Y-, ■ <a u +> Y r ■ <a,,+ ) 


H 


H 


K A .y v m A i3t 

y 4 * ^ ^ 


m Al ,+ 


K K -Yj 

2 4- 1 L -f 


y 2’V 


^ ^ aj |~t ) 


k 3 • r 3 ; 1< 4 - y 3 

+ 3 + 4 

Cm r _ »\^ / - • \ "* 


>'o- < V ) 


1 '4- <a H +) 


*A1 ,+ 


K, 


K, 


K, 


K , 


4 


'3 ^ 2 .a H + r r <V> ? 'o- <a H +> 3 V <a H + -l 
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Som® simplifications can be introduced in the above 
equation. For example, i) the neutral species A1C0H>° can be 
neglected with respect to other charged species. This is mainly 
because A1(0H) J i 3 precipitated as a solid phase under pH 
conditions where A1<0H>° is the dominant species in solution; ii) 
the activity coefficients of the univalent species A1(0H>* and 
AHQH)^ can be considered to be equal i.e. - y 

Therefore , 

m Al j 

3 a Al 3+ / f 

or a A1 ? + = m A / f 


l Al 3+ [ 


K, 


K. 


’W 


r 1 - (a H +> ‘ 


J 'l‘ Ca H + 


) 4 - 


where ' f ’ represents the terms within bracket. 

The values of Kj , K £ and l< 4 calculated from the free 

— 5 004 ** 9 yea 

energy data listed in the Table 6.18 are 10 . 10 ‘ and 

1Q -25 . 166 respect j ve 1 y at 25°C. Using the activity co-efficient 
(y } values from Table 5.6, the function f 1 can be evaluated for a 
given pH to set the relationship between the analytical aluminium 
concentration and activity of Al 3+ ion as suggested by Neal 
(1988). At river water pH value, the dominant dissolved aluminium 


species is Al(OH) 


practically equal to K 4 .Ai <a n +> 


Therefore, the value of ' f becomes 
4 
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Table 6.18 Free energy data CkJ/'mole} for aluminium speciation 
calculation at 25°C and 1 atmosphere 


Reference 


Gibbsite 

: -1154.89 

Robie et al . (1978) 

cahoh> 5 3 



Kao 1 i n l te 

: -5799.56 

—do — 

C A 1 2 S i ^0 ^ 

<QH) 1 

4 


Poorly crystalline : ->766.18 

Kittrick, 1966 

kaol ini te 

CA1 2 Si 2 0«.(0H> 4 ] 


H 4 S|0 4 

: -1507.74 

This report^* 

H 2° 

: - 257.14 

Rob ie et al. (1978) 

Al 5 * 

: - 484.88 

Wagman et al.(1968). 

2 + 

A1 <0H> 4 

: - 695.46 

-do- 

-do- 



A1 (0H>2 

: - 905.50 

Reesman et al.(1969) 

A 1 ( OH ) ~ 

: -1501.25 

-do- 


a Calculated from 6 ppm SiQ^ for quartz solubility 
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The solubility of kaolinite can be expressed 


as 


I A1 2 Sl 2°5 (OH) 4 + 5H+ = A1 ^ + + H 4 Si0 4 + j H £ C 


The equilibrium constant 


K 


kao 1 . 


a Al ,+ ' a H SiO 
4 4 

<a H +> 


lo 9 K kaol. = 


log a M 5. * log a + 3 P H 

4 4 


= log m ftl - log f + log a H SiQ + 3pH 
T 4 4 


log m AlT ♦ log *h 4 S,0 4 = 109 f ‘ 5 pH + 109 K kaol. 


The value of K kaoi calculated from the free energy data listed in 
Table 6.18 can be expressed as log K ka0 ^ = 2.017 at 25°C. The 
above relationship results in a straight line with 1:1 slope when 


[log m 

T 

+ 109 a H . S i 0, 

4 4 

1 ■» 

plotted 

aga i nst 

(log 

f 

- 3pH> in 

F 1 g . 6 . 1 1 . 









Using the 

free 

energy 

data of 

poor 1 y 

crystal 1 i ne 

kao lintte 

listed in 

Table 

6.18, 

a parallel 

straight line 


representing its solubility has the equation 


ri__ m + 1 na a 1 = (log f “ 3pH) + 4.924 

i_l°g m + log a u s jQ J 


4 4 
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Table 6.19 Data for kaolinite stability diagram 


Stn. 

-log ™ a1t 

1 og 

1 . 

6.68 

5.81 

2. 

6.70 

VI 

Ni 

3. 

6.96 

7.87 

4. 

6.92 

7.92 

5. 

6.66 

7.95 

6. 

6.85 

7.14 

7. 

6.70 

6.26 

8. 

6.60 

7.77 

9. 

6.60 

6.18 

10. 

6.77 

6.42 

11 . 

7.00 

8.18 

12. 

6.49 

7.42 

13. 

6.05 

6.74 


' lo9 a H <!in PH 
4 4 


3.63 

7.23 

3.62 

7.21 

3.63 

7.75 

3.50 

7.76 

3.56 

7.77 

3.63 

7.57 

3.61 

7.35 

5.56 

7.73 

5.54 

7.33 

5.58 

7.39 

5.61 

7.83 

5.52 

7.64 

5.62 

7.47 


(log f-3pH) log 


-15.88 
-15.90 
-15.38 
-15.36 
-15.36 
-15.57 
-15.79 
-15'. 42 
-15.81 
-15.75 
-15.31 
-15.50 
-15.67 


109 «h 4 s.oJ 


-10.31 
-10.32 
-10.59 
-10.42 
- 10 . 22 ' 
-10.48 
-10.31 
-10.16 
-10.14 
-10.35 
-10.61 
- 10.01 
- 9.67 


Log m A| + Log 0^5^ 
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-17 -16 -15 -14 -13 -12 -11 

( Log f - 3 pH ) 

6.11 Diagram for kaolinite stability at 25°C in relation t 
dissolved aluminium in water 
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The river water analyses have been utilized to calculate 
the relevant parameters of this plot (Table 6.19). As seen from 
the Fig 6.11, the composition of the Mahanadi river water plots 
along the solubility line for poorly crystalline kaolinite. It 
indicates that this phase may have some influence on the aluminium 
concentration in river water. In fact. Hem (1991). has suggested 
that at certain silica levels, poorly crystallized clay mineral 
species can precipitate from natural waters. 

6*7 Heavy Metals in Sediments and Soils 

Stream sediments and soils contain a host of heavy 
metals derived from the weathering of rocks and anthropogenic 
sources such as mine waste. industrial and domestic effluents. 
Although sediment analyses do not furnish quantitative information 
on the absolute degree of pollution, they can play a key role in 
ascertaining relative factors of enrichment which can trace the 
source of pollution of the aquatic environment (Forstner and 
Wittmann, 1981). A systematic study of heavy metals would, 
therefore, include determination of degree of enrichment with 
respect to some background concentration values. The present 
investigation concentrated on the analysis of less than 45 micron 
size fraction of soils and flood plain deposits to detect a 
general trend. if any. A total of seventeen samples were 
collected from different locations of upper part of the Mahanadi 
basin (Fig. 6.12) and were analyzed for Pb, Zn, Cu and Mn by AAS . 
The method of sample preparation and procedure for analysis are 



161 


given in Chapter 2. 

6. 7. 1 Degree of Enrichment in Limestone Terrain 

Analytical data on Pb, Zn , Cu and Mn in both flood plain 
deposits and soils in the carbonate rock dominated terrain is 

presented in Table 6.20. The relative proportion of these 

elements are generally in the order of Pb < Zn < Cu < Mn . The 
same table also lists the average background heavy metal 
concentration in limestone from different sources for the sake of 
comparison. The enrichment, if any, of these metals due to 

weathering or anthropogenic sources can be best evaluated by 
comparing with background levels in local limestones. Geological 
Survey of India (1989) records a concentration of less than 10-40 
ppm of zinc in limestones of this area while the sediments have a 
range of 59 to greater than 110 ppm with an average value of 57. 
Therefore, it seems that there has been an enrichment of zinc 
which can possibly be ascribed to urban effluents. Samples IF and 
8S with highest values are from Nawapara and Raipur towns 
respectively. On the other hand, lead in the sediments and soils 
(5-41 ppm), seems to be within the concentration range for 

limestone <10-100 ppm). Average copper level in soils and 

sediments is 100 ppm though occasional high values up to 185 were 
also obtained. But this does not seem unusual because Geological 
Survey of India reports a wide range of Cu concentration (<10-200 
ppm) in the local limestone. With regard to manganese, the 
average value of 782 ppm in sediments and soils is less than that 
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in limestones quoted by Drever (1988) from the original data of 
Turekian and Wedepohl. This comparison indicates that there has 
been a relative depletion of Mn. Wedepohl (1978) suggested that Mn 
is more abundant in the carbonate than in the detrital fraction of 
the carbonate rocks. Therefore, weathering of local limestones 
may have caused the depletion of Mn in the residual soils. 

6# 7* 2 Dogr ot Enr l chment in Granitic Terrel n 

Another set of flood plain deposits were sampled in the 
granitic terrain near the source of the Mahanadi river. These 
samples were analyzed for Pb , Zn . Cu and Mn . Table 6.21 lists 
their concentration levels along with the background values in 
local granites and other reference granites from different 
sources. It is observed that average metal values of the 
sediments is in the order Pb < Cu < Mn < Zn . Lead concentrations 
are fairly constant for all the samples analyzed and is in the 
vicinity of the lowest value (20 ppm) reported for the local 
granite. Manganese concentration in the sediments shows 
considerable fluctuation (240-1810 ppm) while the background value 
of manganese in granites listed by Drever (1988) and Krauskopf 
(1979) ranges between 400-450 ppm. A high value of 1810 ppm of Mn 
at Gangrel (sample 18F) situated upstream of the Gangrel dam site 
may be due to precipitation of Mn in the reservoir sediments. 
Copper and zinc are enriched in the sediments. The reason for 
this enrichment is not clear from the limited data available at 
present. However, anthropogenic sources cannot be ruled out. 
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Table 6.21 Concentrat ion(mg/kg> of heavy metals in rocks and flood plain 



depos i t s 

in the 

gran i t i c 

terrain . 



Samp 1 e 







(Flood plain 

Pb 

Zn 

Cu 

Mn 

Relative proportion 

depos i t > 







1 IF 

20 

>1000 

130 

380 

Pb < Cu < Mn 

< Zn 

12F 

31 

60 

125 
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< Mn 

13F 

30 

830 

120 

240 

Pb < Cu < Mn 

< Zn 

16F 

29 

>1000 

125 

260 

-do- 


17F 

<3 
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24 
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< Mn 
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23.1 

-- 

133.5 

646.7 

Pb < Cu < Mn 

< Zn 

Gran i tea 

20-73 

10-73 

10-75 

-- 



Gran i t@# 

17 

50 

20 

450 



Grani te$ 

20 

40 

10 

400 




61 Local granite (Geological Survey of India, 1989) 

# Complied data listed by Drever(1988> 

* Krauskopf <1979) 
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6* 7. 3 Summary 

The overall picture that emerges from the distribution 
and concentration of heavy metals in sediments and soils in the 
upper part of the Mahanadi river basin is as follows. 

a> The relative concentration of heavy metals is Pb < Zn < Cu < 
Mn in carbonate terrain and Pb < Cu < Mn < Zn in granitic terrain. 

b) In the granitic terrain, Pb and Mn values in the sediments 
i and soils are similar to that in bed rock. They could have been 

derived by the weathering of rocks in the catchment area. 
Similarly, in the carbonate terrain, Pb and Cu levels in soil and 
sediment are comparable with bed rocks. On the other hand, Mn is 
perhaps depleted during weathering of the limestone. 

c) Zinc is enriched in samples from both the terrains, more so 
in the granitic region. Urban effluents may be a factor in this 

pattern . 

6.8 Factor Analysis 
6. 8. 1 I ntr odtic 1 1 on 

As stated by Davis (197?). "factor analysis attempts to 
reveal a single underlying structure that is, presumed to exist 
within a set of multivariate observations." This structure can be 
expressed in terms of variances and co-variances between the 
variables and the similarities between the observations 


There 
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are several types or modes of factor analysis possible. R-mode 
factor analysis describes the relations among the variables in 
terms of the samples while Q-mode studies the correlations among 
the sampling sites. Separate interpretations are required for 
these two modes because the final results of such analyses produce 
a set of factors demonstrating the relations. But explanations to 
these relations are to be found from the knowledge of the 
geochemi ca I env i ronment . 

Factor analysis is carried out in 3 steps : 

i) Preparation of the correlation matrix which is produced by 
multiplying a data matrix or a normalized data matrix by its 

transpose . 


ii> Extraction of initial factors so as to explore the 
possibility of data reduction. This is done by diagonalizing the 
correlation matrix and obtaining eigen values and principal 
components (eigen vectors). The factor-1 would then relate to the 
highest eigen value explaining the greatest percentage of variance 
in the data set, factor-2 would explain the greatest of the 
remaining variance and so on. 


iii) Rotation of the initial factor matrix to a final solution 


which provides a set of simple and meaningful factors. 


D i f f erent 


rotations of orthogonal axes derived from principal component 
solution may result in different terminal solutions. In the 
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'varimax* rotation, orthogonal factor axes are located such that 
the variance of each factor is maximized. This results in 
equitable distribution of factor loading. In the 'oblique* 
rotation, the reference factor axes are not constrained to be 
orthogonal. Oblique solutions are appropriate where underlying 
causal factors are suspected to be correlated among themselves. 

6. 8> 2 Fact o* Analysis of Chemical Data of the Mahanadi River 

Wat er 

Water quality data at the 13 CWC stations along the 
Mahanadi river were subjected to factor analysis as discussed 
below. Chemical parameters chosen as variables for R-mode factor 
analysis were SiO,,, Na , K, Ca, Mg, HCO^, Cl and SO^. Monthly data 
for these variables at each sampling site were averaged out on 
annual basis. A set of five analyses for each station resulted 
from the five year data base. The raw data thus obtained were 
normalized by dividing the concentration of each chemical 
parameter by the highest among them (Table 6.22). Normalization 
was felt necessary because i > the basic data lacked normal 
distribution and 1 i )’ variables of higher magnitudes would have 
dominated over those of lower magnitude. Hitchon et al . (1971) 
and Stallard and Edmond (1983) preferred use of logarithmic 
transforms of data to raw data for factor analysis. 

Factor analysis was carried out using SPSS (Statistical 
Package for Social Sciences) with an HP-9000/850 series 
supermini-computer at I.I.T. Kanpur. A minimum eigen value of 0.5 
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6.22 Hormalized data for 


Si °2 

Na 

K 

Ca 

0 

713 

0 0 74 

0 539 

0 312 

0 

536 

0 090 

0 223 

0 264 

0 

554 

0 188 

0 409 

0 363 

0 

431 

0,104 

0.281 

0.270 

0 

813 

0 184 

0 e 46 

0.291 

0 

67 ! 

0 1 06 

0.363 

0.451 

0 

681 

0 173 

0 360 

0.463 

0 

541 

0 171 

0.316 

0.493 

0 

.48 6 

0 . 187 

0,264 

0.387 

0 

166 

0 149 

0.320 

0.417 

0 

677 

0.350 

1 .000 

0.762 

0 

602 

0 483 

0 942 

0.727 

0 

580 

0 , 392 

0 743 

0.774 

0 

440 

0 347 

0 6 9 1 

0 714 

0 

8 36 

0 376 

0 849 

0.759 

0 

916 

0 825 

0 740 

0.935 

1 

0 00 

1 000 

0.909 

1 .000 

0 

592 

0 ©96 

0 635 

0 974 

0 

717 

0 650 

0 498 

0.839 

0 

882 

0 949 

0 648 

0.988 

0 

698 

0 434 

0 749 

0 719 

0 

745 

0 44 7 

0 638 

0 736 

0 

089 

O 463 

0 637 

0 685 

0 

6 64 

0 423 

0 613 

0 719 

0 

861 

OSH 

0 632 

0 6 68 

0 

433 

0 166 

0 287 

0,597 

0 

6 72 

0 270 

0 317 

0,633 

0 

542 

0 283 

0 363 

0 726 

0 

598 
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0 294 

0 587 

0 
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0 351 

0 311 

0 632 

0 
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0 193 

0.430 

0 , 399 

0 

,601 

0,186 

0.415 

0.409 

0 

.593 
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0.399 

0 361 

0 
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0.523 

0 414 

0 

315 

0.317 

0.691 
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0 

837 
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0 

70 S 
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0 
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0 
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0 31 1 
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0 
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0 
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0.41 1 

0,300 

0 

756 

0 113 

0.377 

0.268 

© 

. 68© 

0.110 

0.368 

0.869 

© 

537 

0.109 

0.374 

0.269 

© 

867 

0 117 

0 356 

0,272 

© 

641 

© 13 © 

0 283 

0 328 

© 

78 f 

0 239 

0.233 

0.313 

© . 

638 

0.199 

0.244 

0.329 

© 

60 © 

0 199 

0.840 

0.326 

0 , 

833 

0.847 

0.844 

0.340 

0 

577 

0.327 

0.297 

0.665 

0 

081 

0.476 

0.311 

0.643 

0 

575 

0.485 

0.314 

0.615 

0. 

SIS 

0.407 

0.865 

0.627 

0 

264 

0 526 

0.292 

0.596 

0. 

7 6S 

0 241 

0.368 

0.530 

0. 

790 

0.341 

0.398 

0.478 

0 

704 

0.266 

0.395 

0.498 

0 

688 

0.318 

0,360 

0.494 

0 

859 

0 291 

0.343 

0.480 

0 

471 

0.104 

0.315 

0.475 

0 

571 

0 163 

0 287 

0510 

0 

553 

0 166 

0 299 

0,479 

0 

637 

0 176 

0 263 

0.479 

0 . 

335 

0 176 

0 280 

0,483 


factor analysis Cn=65) 


Mg HCO ? Cl S0 4 

0 166 0 e63 0 . 429 0 . 0 6 t 
0106 0103 0 485 0.125 
0.1 33 0 .204 0.715 0 . 287 
0.181 0. 177 0.400 0.116 
0.126 0.199 0.361 0.153 
0.229 0.356 0.454 0.104 
0.198 0 .302 0.550 0.029 
0161 0.319 0.609 0. 073 
0.146 0.850 0.450 0.116 
0.136 0.270 0.463 0.149 
0.595 0.644 0.804 0.260 
0.516 0 .567 0.740 0 . 254 
0.463 0,577 0.741 0.262 
0.363 0.524 0.637 0.376 
0.449 0,556 0.636 0.897 
0 955 0.955 0.718 0.597 
0.819 0.689 0,750 0.919 
0.720 0.932 0.666 1 .000 
0,675 0.666 0.650 0.563 
1.000 1 .000 0.800 0.621 
0.653 0.626 0.729 0.397 
0 605 0.539 0.900 0.492 
0.643 0 615 1:000 0.468 
0 433 0,534 0.760 0.441 
0.602 0.580 0.762 0.442 
0.329 0.447 0.443 0.157 
0.306 0.465 0.660 0.124 
0,332 0.498 0.622 0.130 
0.886 0 421 0.500 0.149 
0.314 0.473 0.590 0.314 
0.256 0.350 0.517 0,230 
0. IBS 0.283 0.450 0. 179 
0. 161 0.250 0.565 0.899 
0.223 0.278 0.538 0.402 
0.219 0.304 0.612 0.361 
0.436 0.522 0.600 0.274 
0.303 0.373 0.633 0.843 
0.331 0.438 0.716 0.231 
0.339 0.410 0.645 0.811 
0.400 0.469 0.600 0.351 
0.218 0.293 0.436 0.809 
0.161 0.814 0.481 0.171 
0.183 0.213 0.395 0.147 
0. 140 0. 195 0.337 0.303 
0.192 0.197 0.387 0.226 
0.845 0.344 0.421 0.114 
0.176 0.248 0.504 0.113 
0.177 0 . 248 0 . 425 0.151 
0.185 0.260 0.410 0. 357 
0. 162 0.266 0.416 0.196 
0.402 0.574 0.504 0.181 
0 . 377 0 . 548 0.608 0 . 047 
0.354 0.. 575 0.605 0.1 05 
0.364 0.546 0.563 0.140 
0.362 0.543 0.631 0.813 
0,349 0 . 477 0.510 0.166 
0.289 0 . 376 0.436 0.126 
0 . 278 0.406 0 . 554 0 . 053 
0.305 0.420 0.533 0.180 
0.270 0.391 0 . 544 0.201 
0.226 0.397 0.537 0.805 
0 216 0 . 337 0 . 520 0 . 050 
0.210 0.360 0.571 0.006 
0.253 0.341 0.536 0.233 
0 254 0.353 0.533 0,197 
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was assigned for extraction of factors. The size of the eigen 
value represents the variance of the original data that has been 
extracted on to each factor. Principal component analysis 
extracted 5 factors. Commonality was found to be high except for 
Cl<0-85>. and S0^< 0.71). The generally high commonalities 
indicated that most of variance (90%) i s explained by a set of ? 
factors. Initial factor matrix thus obtained by principal 
component solution requires some form of rotation. Oblique 
rotation was preferred to others because the controlling processes 
in hydrogeochemistry are often correlated. Sampson (1968) 
illustrated the use of oblique solutions in the analysis of 
formation water data. Usunoff and Guzman’ (1989) also favoured 
oblique rotation because geochemical processes often evolve 
through very dynamic and highly interrelated mechanisms. 
Therefore, upon oblimin rotation of the initial factor matrix with 
Kaiser normalization, the following factor matrix given in Table 
6.27a was obtained. 


6.8.3 Geochemical Interpretation 

The factor analysis shows three trends or groupings 

among the variables. They are : 


i ) 

Na-Ca-Mg-HCO ? 

: 

Factor-1 

i i ) 

Si0 2 

: 

Factor-2 

i i i > 

K-Cl 

: 

Factor-? 
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Tabl© 6.23a Obi i mi n rotated factor matrix of chemical data 

C R-mode) 


Var i abl a 


Factor 


Communal i ty 


1 

2 

3 


sio 2 

0.160 

0.996 

0.149 

0.993 

Na 

0.968 

0.120 

0.629 

0.945 

K 

0.651 

0.178 

0.946 

0.900 

Ca 

0 . 9 ?8 

0.1?? 

0.768 

0.901 

Mg 

0.967 

0.221 

0.745 

0 . 949 

HC0 ? 

0.971 

0.191 

0.679 

0.946 

Cl 

0.744 

0.042 

0.909 

0.852 

SCL 

4 

0.837 

0.241 

0.624 

0.710 

Eigen value 

5.616 

0.986 

0.594 


(Principal 
component 
solut i on) 





Percentage 
var i ance 
expla i ned 

70.20 

12.30 

7.40 


Cumulat i ve 
percentage 
of variance 

70.20 

82.50 

89.99 



Table 6.23b Factor correlation matrix 


Factor-1 

Factor-1 1.00 

Factor-2 0.1? 

Factor-? 0.71 


Factor-2 Factor-? 

1.00 

0 . 1 ? 1-00 
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Factor-1 : More than 70% of variance among the ions is 
explained by this factor. The interrelationship between Ca, Mg 
and HCO ? in water can be easily explained by carbonate weathering 
in the basin. On the other hand, association of Na with these 
ions may be related to the weathering of silicate impurities in 
limestone. As discussed earlier, the acid-leached residue of the 
do 1 om i t i c 1 i me stone contains albitic felspar. Therefore, factor — 1 
can be interpreted as essentially carbonate weathering factor. 

Factor-2 : This factor explains more than 12% of total 
variance among the ions and is characterized by high factor score 
for silica. Therefore , this factor may be interpreted in terms of 

weathering of silicate rocks. 


Factor-5 : This factor accounts for only 7% of total variance 
among the variables. Significant loadings were obtained for K and 
Cl followed by moderate loadings for other ions. Although a 
control of ©vapor i te weathering is indicated, available geologic 
information is insufficient to confirm this possibility. This 
factor, however, is related to factor-1 (r = 0.71) as obtained 
from the factor correlation matrix (Table 6.25b). 

Another R-mode factor analysis with the same variables 
was also performed. Initial factor matrix was subjected to 
* var i max * rotation with Kaiser normalization. The varimax 
converged in four iterations. The rotated factor scores for 5 
factors extracted (Table 6.24) showed no significant changes from 
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Table 6*24 

Varimax 

rotated 

factor matrix CR-modeD 


Var i abl e 



Factor 



i 

2 

3 

SiO £ 


0.079 

0.059 

0.992 

Na 


0.933 

0.267 

0.057 

K 


0.333 

0.872 

0.121 

Ca 


0.823 

0.469 

0.067 

Mg 


0.866 

0.417 

0.155 

HCO- 

> 


0.907 

0.329 

0.127 

Cl 


0.302 

0.775 

-0.019 

S0 4 


0.769 

0.342 

-0.034 

Percentage 




■ 

variance explained 

70.20 

12.30 

7.40 

Cumulative 

per- 




centage of 

variance 

70.20 

82.50 

89.90 
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the obi iquely rotated factor scores. But factor-2 and factor-5 in 
the oblique transformations were interchanged i.e. factor-2 is now 

related to K-Cl while factor-5 is Si0 2 factor. The percentage of 

variance explained by K-Cl is 12 and that of SiO is 7. This 
indicates that the basic composition of factors did not change. 
Therefore, the variance of chemical constituents in adequately 
explained by orthogonal and oblique factors but the latter was 
preferred for geochemi cal interpretations. 

The primary control of factor-1 (carbonate weathering) 
on Mahanadi river water composition is clearly indicated by the 
statistical method adopted in this section. This agrees with the 
conclusion drawn from the ternary diagram;? developed in the 
preceding section which were based on standard geochemical 


arguments . 


CHAPTER-7 

WEATHERING IN PILOT AREA 


7.1 Introduction 

Th. ult.rn.te product of rock-water interact. on discussed 
th. previous chapter is an assemblage of secondary minerals 
which appear in the soils and sediments. The formation of soil 
over limestone bed rock is particularly interesting. The 

conventional scheme of limestone weathering involves congruent 
dissolution of carbonate minerals. Theoretically. therefore. a 
pure limestone bed rock cannot support a residual soil profile. 
One would expect to see only solution features and a karst 
topograpny . On the other hand, most limestones contain enough 
silicate minerals to give rise to clay-rich soils which normally 
develop over silicate terrains. The study of these soil profiles 
requires large scale mapping of a small area where all information 
regarding limestone weathering including mineralogi cal and water 
chemical data can be readily obtained. For th 1 s purpose . the 

o » o * 

terrain around Baloda Bazar <21 59 N , 82 1 0 E) was selected as a 
Pilot Area for recording the pattern of soil formation over a 

stromato 1 1 1 1 c limestone bed rock. 
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7.2 Geology of the Area 

The formations exposed in the area belong to Raipur 
series. From the geologic map (Fig. 7.1), it is seen that 
Tarenga shale. Chandi limestone and Gunderdehi shale are exposed 
in the pilot area. Stratigraphic superposition and primary 
minerals of these rocks have been discussed in Chapter 5. 

Limestone rocks belonging to the Chandi formation cover 
a considerable portion of the pilot area. These stromato i i t i c 
limestones are indicative of a shallow marine condition of 
sedimentation (Murti, 1987). The weathering surface of limestone 
shows solution features like open cavities, solution caves, open 
joints, stylolites etc. Because of their deposition in a near 
shore environment, these limestones contain a relatively large 
proportion of silicate minerals in addition to common carbonate 
minerals. The weathering of this impure limestone should involve 
parallel congruent and incongruent reactions. An attempt has been 
made to interpret the mechanisms of the weathering of carbonate 
and silicate fractions and consequent soil formation. 

7.3 Weathering of Carbonate Minerals 

Calcite and dolomite are the major carbonate minerals 
present. In addition, ankerite was detected in some samples on 


the basis of the characteristic X-ray peak between 2.89 and 2.9 A. 
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It wa* conf i rmed by the presence of an endothermic OTA peak below 
800 °C <F,g. 7.2). Weathering of these carbonate minerals takes 
place by a process of congruent dissolution in presence of 

C0 2 -saturated water: 

a) CaCQj + C0 2 + h^O = Ca 2+ + 2HC0~ 

b) CaMg(C0 3 > 2 + 2CQ 2 + 2H 2 0 = Ca 2+ + Mg 2+ + 4HC0‘ 


As the mineral goes completely into solution, the 
limestone develops cavities, sink holes and cavernous zones 
depend s ng on the scale of weathering. If C0 2 escapes from the 
solution saturated with calcite and dolomite, the above reactions 
are reversed and r epr ec i p i tat i on occurs. Cave deposits and 
carbonate concretions are well known examples of this situation. 

A set of groundwater analyses from this area (Ramanna, 
1976) was utilized to test saturation with calcite and dolomite. 


The calculations summarized in Table 7.1 show that the ionic 

-8.82 -7.46 

activity product . a^g 2- range from 10 to 10 and 

“ 3 

brackets the solubility product of calcite, = 10 ‘ at 25°C. 

2 

Similarly, the value of ag a 2 + * a Mg^ + " ^ a C0 3 2- ^ ranging from 


10' 17 - 76 to 10' 


15.23 


brackets the solubility product of dolomite. 


K a 10 - 17 -° at 25°C. Most of the ground water samples as well as 

d 

Seonath river water are super-saturated with respect to calcite 
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T,i! 1 «’ 7.1 Gi 1 c u I «i t ed ion activity products in river and ground 
water of Chhatt i sgarh Basin. Original water analysis 

from Ramanna ( 3976). 


Sample No. 
i 1 oc a 1 1 cn 


leg ta „ .a 7 ] 
Cu COf" 


losla r 2 + ' a 

Ca Mg CO^ 


IC - 3 
i Bert a ra ) 

-7.88 

- 15.50 

1C - IS 
(Charm r 1 3 

-7.46 

- 15.26 

X - 6 
(Dhobnidih / 

-8.44 

- 17.02 

X - 9 
<Kef.li ) 

-8.82 

- 17.57 

•V* - 1 

* 

(Diggi 3 

-7.84 

- 15.38 

X - 3 

(Bhatbera ) 

-8.53 

- 17.76 

X - 12 
(Bltkul l ) 

-7.80 

- 15.57 

X - IS 
( Senodha 3 

- 7.72 

- 15.70 

1C - 5 
(Gogia ) 

- 7.67 

- 15.23 

Seonath River 
(Nandghat ) 

- 7.67 

- 15.28 
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and dolomite. This may explain tho 

t "Hum une widespread occurrence of 

carbonate concretions in this area. 

7.4 Weathering of Silicate Minerals 

Limestone bed rocks are observed to support two types of 
soil cover which match with the calcimorphic soils described by 
Loughnan (1969). One is a concretionary red soil (Terra Rossa) . 
The other is a plastic, organic rich black soil (Rendzina). X-ray 
diffraction studies (Table 7.2) revealed that the assemblage 
kaol mite * illit® + quartz is common to both soils. However, the 
red soil is characterized by concretions of goeth i te while a 
swelling clay ( montmor i 1 1 on i te ) and poor 1 y crystalline chlorite 
were detected in some samples of the black soil. This soil 
mineralogy can be traced to the silicate minerals trapped during 
the deposition of the stromatoi ltic limestone. Examination in 
hand specimen and thin section along with X-ray pattern of HC1 
leached residue of this limestone (Fig. 7.5) showed the presence 
of quartz. alka.i felspar, kaolimte, illite, chlorite and 
glauconite. Among these glauconite is probably a primary source 
of iron for the *'ornation of goethite by weathering as discussed 
in the following sub-section. Additional iron can be derived 
i n the soil zone by dissolution of ankerite. 

(Ca, Mg. Fe)C0 5 + CO., + H.,0 ► Ca 2+ + Mg 2+ + Fe 2+ + HCO^ 



Table 7.2 Soil mineralogy in the pilot area from XRD 
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Red soil (Ter ra Rossa) Black soil (Rendzina) 


■m 


Mineral ogy 

Sample No. 


Mineralogy 

S-5 

Q 

♦ I + K + Go 

S-4 

Q 

+ I + K + Mo 

(Damru > 



(Tarasiv) 



S-6 

Q 

*■ l + K + Go 

S-8 

Q 

+ I + K + Chi 

(Sonbar sa ) 



(Karmada) 




Q-Quartz . I - I 1 1 i te . K-Kaol ini te , Go-Goethi te , Mo-Montmori 1 Ionite, 


Chi -Ch 1 or i to 
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4* 1 Identification and Weathering of Glauconite 

Identification of glauconite in limestone was carried 
out through a series of steps. Examination in hand specimen and 
thin sections showed green coloured, rounded to elliptical 
pellets. Particles in thin section have optical properties 
similar to chlorite. However, the X-ray diffraction pattern of 
acid-leached residue (Fig. 7. 3) showed a sharp 10A°peak which 
could be due to a mica-type mineral. Therefore, in order to 
confirm the presence of glauconite, such green coloured pellets 
were separated from the carbonate matrix for analysis under 
scanning electron microscope. SEM photographs showed laths of 
irregular crystals. Energy dispersive X-ray analysis (EDAX) of 
these grains (Table 7.3> established them to be potassium and 
iron rich silicate. The composition provided a fairly good match 
with the idealized glauconite given by Burst (1958). 

The weathering of glauconite has been recorded by 
Loughnan (1969) and Courbe et al . (1981). In the soil 


env i ronment , 

g 1 aucon i t e 

initially alters 

to 

a montmor 1 1 1 on l te 

( nontron i te ) 

which in 

turn gi ves rise 

to 

an assemblage of 


kaolinite and goethite. Wolf (1967) suggested that spheroidal 
particles of goethite could be pseudomorphs after glauconite. It 
is obvious that these mi neral og i cal changes are similar to 
standard silicate weathering reactions with the added parameter of 



Table 7.3 Chemical composition of glauconite by SEM-EDAX 


El ement 

Normal i zed wt . % 

G-l G-2 

Ideal glauconite 

Na 

0.00 

0.00 

1 

\ 1.84 

Ca 

9 .60 

8 . 34 

J 

f 

! 

Mg 

0 . 00 

0.00 


5.67 

A 1 

11.93 

15.40 


13.90 

Si 

33.71 

44.28 


48.51 

K 

19.91 

17.73 


14.18 

Fe 

22.83 

14.25 


19.29 

Total 

100.00 

100.00 


100.00 

G-l , G-2 are 

glaucon i te 

pellets in 

1 imestone . 



♦Element composition recalculated from the glauconite formula 
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oxidation of ferrous to ferric iron. The ferric iron is 
ultimately fixed in goethite. 

From a crystal lo-chemi cal point of view, glauconite is a 
ferrous-ferric equivalent -of illite. An idealized formula adopted 
from Burst (1958) can be written as K, _ ( Fe^ + , , Fe 2+ „ , >X where X = 
A1 g 66 Si 3 74°20 <OH) 4' Weathering in the presence of dissolved CO^ 
is essentially a reaction with H + ions contributed by the weak 
acid H^CQ^ • Conversion of glauconite to an ideal nontronite. 
Kg ^^(Fe^ + )^X involves the following steps : 

glauconite = K. ,,F ej* X + 0.84 Fe 2+ + 0.84 K + 

0.665.16 

0.84 Fe 2+ + 0.21 D £ + 0.84 H + = 0.84 Fe 5+ + 0.42 H £ 0 

K n ..Fe^ + W X + 0.84 Fe^ + = nontronite 

0.66 >.16 

glauconite + 0.21 + 0.84 H + = nontronite + Q.84K* + 0.42 H^G 

Similarly, the weathering of nontronite to kaol mite 

plus goethite involves release of K + ion- and dissolved silica, 

H SiQ . The overall reactions for giauconi te-nontronite, 

4 4 

non t r on 1 1 e -kao 1 i n i t e + goethite and glaucon i te-kao 1 i n i te +goethite 
euilibria have been provided in Table 7.4. The Gibbs energy 

values used to calculate equilibrium constants of these reactions 
have been listed in Table 7.5. From these data, a stability 




Table 7.4 Equilibrium constants for stability of glauconite in 
the soil environment at 25°C, 1 atm. 


Glauconite (gl), ^ ^(CH^ 

Nontronite (no). K 0 ^(Fe 3 *)^ ggSi, . 34 O 20 (OH^ 

Kaolinite (ka), A1 , Si ,0. ^ (OH ) Q 

4 4 JLU o 

Goethite (go), FeOOH 
gl-no boundary : 

gl + 0.21 0 2 + 0.84 H + = no + 0.84 K + + 0.42 H 2 0 

log (a /a ) 

K H 

= 19.3 for P n = 0.21 atm in earth's atmosphere 

2 

no-ka+go boundary : 

6 no + 4 H + + 82 H 2 0 = ka + 24 go + 4 K + + 40 H 4 SiC> 4 

l°S<a /* > + 10 log a H,SiO, * ~ l3 - S5 


gl-ka+go boundary : 

6 gl + 1.26 0 2 + 9.04 H + + 79.48 H 2 0 

= ka + 24 go + 9.04 K + + 40 H 4 Si0 4 

log ( a /a ) + 4.42 log a - 

tv H 4 4 


-8.64 


Note - Relevant Gibbs energy data are listed in Table 3. 
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diagram (Fig. 7.4> in terms of the thermodynamic activity of K + , 
H + and H^SiO^, has been constructed by the method of Garrels and 
Christ (1965). The idealized composition might introduce 
considerable uncertainty in the position of boundaries. However, 
from the position of the glaucon i te-nontron 1 te-kao 1 1 n i te triple 
point, it is reasonable to conclude that glauconite weathers to 
kaolimte plus goethite at relatively high pH, high K + ion 
activity and low silica activity. 

7.5 Relative Mobility of Elements 

The proposed weathering scheme for the carbonate and 
silicate fractions in lime stone is based on the removal of major 
elements in solution and fixation of A1 in clay minerals and Fe in 
goethite. Comparison of chemical analysis of limestone bed rock 
near Baloda Bazar and of the local soil confirms the same trend. 
Table 7.6 shows the oxide analysis recalculated to element 
percentage. The soil to rock ratio is high for Ai and Fe but ■ much 
lower for Ca and Mg. These values indicate that the relative 
mobility of elements during weathering has been Ca > Mg > Na > K > 


Si > Al >Fe . 
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Table 7.5 Standard state Gibbs energy of format ion of minerals 

and related species for glauconite weathering reactions 


Substance 

AG°f 

KJ /mole 

Sou rce 

Glauconite 

-9296.80 

1 

Nontronite 

-9053.30 

1 

Kaolinite, 

- 7598.72 

2 

A1 4 Si Ao (OH, S 



Goethite 

-488.55 

2 

° 2 S as 

0.0 

2 

H + 

' 0.0 

2 

K + 

-282.49 

2 

h 2 o 

- 237.13 

2 

H.SiO. 

4 4 

-1307.74 

3 


(1) Calculated by the oxide summation method (Tardy and Garrels, 
1974). 

(2) Converted to KJ/mole from Robie et al. (1978). 

(3) Calculated from 6 ppm Si0 2 for quartz solubility and AG°f 

of quartz = -856.29 KJ/mole (Robie et al., 1978). 



tog a n 4 Si0 4 

Stability fields of glauconite (gl). nontronite (no 
and kaolinite (ka) with goethite (go) at 25^- 
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Table 7.6 Calculation of relative mobility of elements during 
weathering of limestone in Bhatapara area 


A. Chemical analysis (TISCO/Sonadih Project, unpublished report) 


Average limestone Average soil 


Si0 2 

7.72 

39.52 

A1J0 

2.91 

16.57 

F '2°3 

0.98 

27.20 

CaO 

46.89 

1.78 

MgO 

2.40 

0.72 

NaJD 

0.14 

0.35 

K 2° 

0.33 

0.82 

LOI 

38.24 

13.69 


99.61 

100.65 

Element concentration 

(Recalculated to 

100 percent ) 



Limestone 

Soil 

Soil /Limestone 

Si 

8.77 

37.80 

4.31 

A1 

3.73 

17.86 

4.79 

Fe 

1.66 

38.90 

23.43 

Ca 

81.40 

2.60 

0.03 

Mg 

3-51 

0.89 

0.25 

Na 

0.25 

0.53 

2.12 

K 

0.66 

1.39 

2.11 




CHAPTER 8 


CONCLUDING REMARKS 


8. 1 Summary 

This thesis describes the geochemical aspects of 
rock weathering in the Mahanadi river basin. A discussion on 
weathering in major world river basins and possible controls on 
river water composition is followed by a review of Indian river 
basins. This leads to the objectives of the present work. A 
representative area of bed rock exposures in the upper reaches of 
the Mahanadi river has been selected for detailed observation. 
The effect of recent deposits and human activities is minimum in 
this area. Moreover, the exposures of silicate and carbonate 
rocks were expected to provide a contrasting picture of rock 
weather i ng . 

Geomorpho 1 og i cal analysis of a tributary in the source 
area indicates that the river is actively eroding granitic rocks. 
The hypsometric curve shows that 64% of ' the land mass has been 
eroded out. On the other hand, a study of satellite imagery and 
field observations reveals the presence of karst features in the 
adjoining limestone terrain of the Chhattisgarh basin. 

Size analysis of suspended sediments of four upland 
tributaries show that 70% of the material is carried in the silt 
fraction. These sediments contain unweathered felspar and 
occasional calcite. These features point to a transport dominated 
(weathering limited) denudation regime in the upper part of the 
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river basin. Sediment load is strongly related to discharge 
because calculations show that nearly 99 % of the suspended load is 
carried during the monsoon period. Total annual load of the 
Mahanadi river at the farthest downstream CWC sampling point near 
Tikarpara is 35.24 million tonnes out of which physical load 
constitutes 30.89 million tonnes and chemical load is 4.35 million 
tonnes. The overall ratio of physical to chemical erosion up to 

this part of the basin is estimated to be around 7:1. 

2^ 4" 2 -t* -f* 4* — — 

Ca , Mg , Na , K are major cations and HCO, , Cl and 

2 - 2+ 2+ 

SO, are major anions in river water. Ca , Mg and HCO, show 

4 3 

lowest concentrations during the monsoon period. Dissolved silica 
and sulphate concentrations are directly related to discharge. 
The concentration of major ions at the confluence of tributaries 
matches fairly well with the predicted value obtained by 
area-weighting the concentration in individual tributaries. This 
suggests a substantial contribution from weathering of rocks in 


the catchment area. 


Dissolved CO, in river and ground water acts as the 
principal weathering agent. Apparent partial pressure of * n 
the river water is more than the atmospheric CO, pressure. The 
acquisition and proportion of three major ions Ca, Mg and HCO^ in 
river water have been evaluated in terms of weathering of 
carbonate and silicate minerals in the rocks. Using the global 
average for fractionation of Ca and Mg (Holland, 1978), the total 
HC0 5 in river water can be divided between parts derived by 
carbonate and silicate weathering. It is found that the ratio of 
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carbonate derived HCO ? to silicate derived HCO ? , approaches a 

value of 2.0 for the limestone terrain while it is more than 1.0 

in other parts of the basin. This implies a relatively greater 

influence of carbonate weathering over silicate weathering even in 

the downstream part of the river basin. This conclusion is 

confirmed by a balance between (Ca+Mg)^ . , and (HCCLK . , as 

total 3 total 

well as a cluster of river water analyses near the HCO^ vertex in 

a ternary plot of HCQ,, SiO^ and SQ^ on an equivalent basis. 

Further support comes from an R-mode factor analysis of major ions 

in water in which Na-Ca-Mg-HCO^ factor alone explains 70^ of the 

total variance among the ions. 

Carbonate and silicate mineral equilibria have also been 

evaluated in terms of thermodynamic data. The ion activity 

_ s 194 

products for calcite and dolomite fall in the range of 10 ' to 

10 9.748 and 16.494 to jq 19.918 respec t i ve 1 y , indicating a 

general undersaturation with respect to these carbonate minerals in 

river water. This explains the aggressive weathering of dolomitic 

limestone in the river basin. Water analyses from different 

sampling stations plot within the kaolinite field in the standard 

stability diagrams for K, Na and Mg systems. On the other hand, 

the analyses plot near the kaolinite-Ca montmor i 1 Ion i te boundary 

in the stability diagram for Ca-system. This suggests that 

kaolinite and Ca-montmor i 1 1 on i te can co-exist in the river water 

environment. This theoretical observation is supported by the 

o 

presence of kaolinite and a glycol expansive 14A clay mineral in 
suspended sediments and soils. The dissolved A1 concentration 


195 




also indicates the stability of poorly crystalline kaolinita phase 
in the river water environment. 

Analysis of sale 'ted heavy metals in soils and sediments 
show a general distribution pattern of Pb < Cu < Mn < Zn in the 

granitic terrain. Zinc specifically, seems to be enriched in 

. 

soils and sediments from both the terrains indicating a possible 
anthropogenic source. 

A detailed mi neral og i cal study of dolomitic limestone 
interbedded with shales and the pattern of soil formation in a 
pilot area within the Chhattisgarh Basin indicate two parallel 
weathering schemes. A congruent weathering of calcite and 
dolomite results in typical solution features and incongruent 
weathering of glauconite gives rise to a red lateritic soil cover 
over the limestone bed rock. The order of relative mobility 
during weathering of elements is Ca > Mg > Na > K > S i > A1 > Fe. 

8-2 Suggestions for Further Work 


1 > 


As discussed above, the thesis project focussed on the 
weathering of bed rocks exposed in the upper part of the 
Mahanadi basin. There were indications that heavy metals 
like Zn in soils and sediments may have an anthropogenic 


source. This aspect can be further elaborated. A correlation 
with mining areas and high Fe, Mn , Cu, Zn and Pb 
concentration has been suggested by Chakrapani and 


Subramanian (1990). 


The Hasdo. Mand and lb valleys may be 
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suitable pilot areas to decipher the role of coal mining and 
chemical industries in contributing trace metals to river 
water and sediment. 

2> The discharge and sediment transport by the Mahanadi is 

controlled by one of the largest earth dams in the country 
located at Hirakud. A comparison between upstream and 
downstream segments through controlled sampling should bring 
out interesting patterns. 

5 ) There have been reports on the possible impact of polluted 
river water on the construction materials of the Hirakud dam. 
This aspect is worth future investigation. 

4) The influence of alkaline soils and saline groundwater on 

river chemistry can be studied particularly in the alluvial 
deposits near the river mouth (Sarin et al . , 1989). 

5> The role of rock-water interaction has been interpreted in 

terms of thermodynamic equilibria in this thesis. It would 
be important to go into the kinetic aspects of these 
reactions. An attempt may be made to estimate rates from 
weathering rind, as outlined by Whitehouse et al . (1986). 

6) For a complete mass balance study, the influence of bio-mass 

on uptake of weathering related chemical constituents must be 
documented. A study in small watersheds similar to the works 
of Cleaves et al . (1970) and Velbel (1986) may be taken up 

within the Mahanadi basin. 

7 ) Both stable and radioactive isotopes have been successfully 


uti 1 ized in river 


geochemistry in recent times. 


A 
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18 2 

distribution pattern of 0 and H may supplement the 

weathering models established in this thesis. Sedimentation 

226 210 

rates may be determined by the Ra - Pb method on cores 
recovered from river banks (Subramanian et al . . 1985). 
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CAppendix I> 

C A C T I C A L } 


C THIS PROGRAM IS WRITTEN IN TURBO PASCAL. IT ESTIMATES MOLE/LITRE, IONIC ) 
C STRENGTH, IONIC ACTIVITY COEFFICIENT, ION ACTIVITIES FROM THE GIVEN PPM ) 
C CONCENTRATIONS OF IOMS IN WATER. IT FURTHER CALCULATES ION ACTIVITY } 
C PRODUCTS FOR CALCI TE , DOLOMITE AND GYPSUM. } 

program pan l grahy < l npu t , outpu t ) ; 
conit 

m*x_*xp ■ 1 0 0 1 

max^com a 20 t 

A ■ 0„ * 

fcypa 

name * ■ tr ing[4 ] ; 

anal ys I s^arr * array C 1 . .max^exp , 1 .. max_comJ of real; 
comp^arr * array C 1 . .max_com] of real; 

name_arr * array [ 1 . .max__com] of name; 

var 

pH * arrayCl ..max_exp] of real ; 
ac t_va 1 , ppm i ana 1 y * i s_ar r ; 
valency, fW t cormp_arr; 

lonjumi t nami^ar r t 

( , i , tot_ana . 

tot_comp t integer; 

Ionic a tr , I on i c__co«f f , l og_coe f f , m 


j real? 
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G.H r t,xt« 

procedure p r i n t_jna t r i x_va 1 u e ( matrix t ana 1 y s i s_arr ) ; 
vtr 

i,| i i nt agar t 

bag in 

for I i* 1 to tot_ana do 
bag in 

wr i t a ( 1 , * * > i 

for ) i" 1 to tot_comp do 

wr l t a (mat r l xC l , | 3 t 5 » J , * ’); 

wr l t« 1 n i 

and i 

and i 

procadura read — ion^lnformit ion (max^val t integer i Var i on_name t name_ar r i Var val,fwgt i comp^arr) j 
vmr 

I * Integer! 

begin 

< for l i* 1 to max_val do 

begin 

readln(G. i on^nameC i 3 , val C i 3 , fwgt C i 3 ) i 

and t 

and i 

procedure read^ppm <max_vall, max_vaI2 s integer; var ppm r anal ys i s_arr ) i 

I 
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var 

I « ! i Integer* 

begin 

for l t* 1 to max__va 1 1 do 
beg i n 

wr i te InC * Enter values for Analysis No t ',i>« 
wr i te < * P 1 ease give the pH value i '); 

readln<pHC i 3 > ; 
for J i * 1 to max_va 1 2 do 
beg i n 

write(*Enter ion concentration in ppm ' , i on^nama t i 1 , * i ’>* 
read 1 n ( ppmt i , j 3 ) * 
end ; 

wr i t a 1 n ; wr i te 1 n ; 

end % 

end i 

function mo 1 e_pa r_l i t ( f un_ppm, f un_fw t real) i real* 

begin 

mo 1 e_pe r__l i t : = ( fun_ppm/ fun_fw ) * 1 E - 3 : 

•nd< 

function cal cu 1 at e_i on i c_s t r ength (ana_r»o : integer) t real; 
var 

l t integer; 


value ; rea 1 ; 
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begin 

va I u« y * '0 i 

for I i - 1 to tot_comp do 

value «» value ♦ mol e_per_i I t <ppmCane_no, i 3 , fwC 1 3 > * vtltncyCi] * vtl«ncy[i]i 
Ca l culate_i on 1 c^mt rangth i * 0.5 * value; 

•nd t 

procedure co3_act (bi carb_act ,pHi i real: var co3 t real): 
begin 

co3 t * < b I carb__ac t *exp< < -1 0 . 33 ) *ln(l 0 ) > / < exp( -pH i #ln< 1 0 > > ) > * 

endi 

Procedure calculate_lon_act l v l ty < var act_va I t ana l y s i s__arr ) t 
vfir 

1,1 i Integer* 

l apical , lap^dol , lap_gyp,co) i real ; 

begin 

for I i * l to to t_ana do 
begin 

wr t telnCH. * Anal ye i a No. t ’,1); 

wr 1 1 e 1 n < H , ’ 1 on ppm mole/lit ionic-str Ionic act. coeff 

wr 1 te I n ( M . * 

for | i« 1 to tot_comp do 
beg i n 

lonic^itr ** ca 1 cu 1 at e_J on i c_e t rength ( 1 ) i 


! 


ac 1 1 v l ty * ) t 



I og^coe f f *» -1 * A * valancyCj} * valencytjj * 

< »qrt < I on i c_® tr> / < 1 + ®qrt ( Ion l c__® tr) ) - 0.2 * ( I on i c_s t r > > » 
m i* mo 1 »__per_l i t ( ppmC i , j 3 . f w[ j 3 > ; 
lonl c_co»f f i* EXP( 1 og_coef f * ln(10))i 
ac t_va l [ i , j 3 i* m * I on i c_coef f ; 
wr l t® (H, l on_nam®C } 3 , * *>i 

wr i t®(H,ppmC i , j:jt7r3 % * *> ; 

wr l t®(H*m, * * ) I 

wr i t® <H, I on i c_* t r i 5 * 4 , ' *>; 

wr i t« (H, Ion i c_co«f f i 4 i 3 , * ? >; 

writ® <H t act_val C i , J3 > ? 
wr I ttln(H) i 

• nd t 

i t®ln(H, * - - 

wr l t® In(H) t 

wr I t® in(H . * I . A. P® ar® t- * ) j 

wr l t®in<H, f * > i 

writ®<H, 'Cmlcita* >i 
co3_act(act_valU , 3 3 , pHC i 3 . co3 ) ? 
t ap_ca 1 t* co 3*ac t_va 1 C i .33: 
wr 1 1® InCH. 1 * t iap_cal > i 

wr l t® ln(H) i 
wr ( t • ( H , * Do 1 om I t • # ) < 

i ap_do l i* ®ct_valt l * 33*act_jval[ i ,43*co3*co3; 
wr 1 1 • 1 n ( H • * ’ » l ap_do 1 ) j 

wr i t® InCH) t 


wr 1 1® <H . * Gyp«um * ) t 
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/ 


and; 


iap_gyp x * actual C i • 5 3*act_val C i , 71 ; 
wr i t • 1 n < H , * • , iap_gyp> ? 

wr i taln<H, ’ * > ? 

wr i ta 1 n(H) ? 

and ; 


begi n 

am a * gn< G, ' i on . data ’ > i 
reset (H , ’output*); 

wr i te < * Enter maximum number of analysis' : ’>; 

readl n< tot_ana) ; 
wr i t a l n ? 

wr i te ( * Enter maximum number of ions in each analysis : ’); 

read 1 n ( to t_ comp) ; 
wr i te 1 n ; 

wr i ta In < 'Now, read the input values for ions.*); 

read_ion_ informat ion (tot^comp, ion_nam«,valancy,fw) ; 

wr i ta ln( ’Now, read tha values of ppm for each analysis.’)? 

read__ppm< tot_ana . t ot_comp , ppm) ; 

ca l cu late_i on_act i v i ty < act__va 1 ) ; 

close (G) ; 

close <H> ? 


end . 
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(Table 1) 

C ION DATA } 
{ ******** } 


Na 

1 

22.99 

K 

1 

39.1 

Ca 

2 

40 . 08 

Mg 

2 

24.31 

HCO t 

> 

1 

61.018 

Cl 

1 

35.45 

S0 4 

2 

96.06 

Fe 

3 

55.85 

A1 

5 

26.98 

NH. 

4 

1 

18.01 

NO j 

1 

62.01 

P0 4 

3 

94.97 

SiO ? 

0 

92.00 
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APPENDIX II 

C Summarized from Holland, 19781 

Cll Contribution of Sedimentary Rocks to IDS in river water 

Out of total land surface of the earth, igneous and 
metamorphic rocks = 24.535 and sedimentary rocks = 75.7% 

Also, TDS from sedimentary terrains 

* 2 x TDS from igneous and metamorphic terrains. 

If X is the contribution from igneous and metamorphic rocks, then 
weighted fractions for sedimentary and igneous + metamorphic rocks 
are 75.7 x 2X and 24.5 x X respectively. 

Therefore, the ratio 

sedimenta ry _ _75.7 x_2X = £, 25 

igneous + metamorph i c ~ Z3T:T~x W~ 

For TDS. 

sedimentary + (igneous + metamorphic) = 100 

or, 7.25 (igneous + metamorphic) = 100 

or, (igneous + metamorphic) = = 13-85 = 14% 

Therefore , 

the contribution of sedimentary rocks = 100-14 = 86% 


C2> Ca and Mg in Sedimentary Rocks 
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Define X 


Ca in carbonate minerals 


Ca 


Ca 


tot 


where all concentrations are in moles/kg of rock. 

sulphate minerals, CaO = CaO - SO, 

corr > 


X 


max 

Ca 


CO, 


Ti0‘ 


corr 


when CO^ < CaO 


= 1.0 when CO^ > CaO 

CO. 

t mi n 2 

Ca = CatT" 

corr 


when 




S 


CaO 

CaO 


car 

tot 


b 


+ MgO 
+ MgO 


carb 

tot 


Available data from Russian and American sources 
Tables 4-12) show that for platform sediments. 


max 

Ca 


1 . 0 , 1.0 


,mi n 
Ca 


5E 


0.81, 0.87 


aV9 *Ca 


0.92 ± 


Correcting for 


(Holland, 1978; 


0.08 



?13 


for geosynclinal sediments , 

avg X Ca - 0.80 ± 0.10 

Existing data also indicate that the proportions of mass of 
platform and geosynclinal sediments are 2 7 % and 7 1 % respectively. 
Therefore, weighting the Ca and Mg fractions accordingly. 


X Ca * 0.62, 0,81 


X" = 0.85. 1.0 

Ca 



<0.92 ± 0 . 08 > x 0.27 + <0.82 ± 0.10) x 0.75 
0.85 ± 0.08 




214 


For geoaynclinal sediments. 


avg X u = 0.45 ± 0.20 
Mg 

.*• X Mg = <0.64 ± 0.20) x 0.27 + <0.45 ± 0.20) x 0.75 

* 0.50 + 0.20 


Not# s 1> TDS = Total Dissolved Solids, corr = corrected, 
maximum, min * minimum, avg = average, carb = 
tot * total . 


Mg 


,mi n 
'Mg 


0.62, 0.81 


0.00, 0.56 


max = 
carbonate . 


2) Compared to Ca, more Mg is tied up with silicates and 
uncertainty of Mg value is greater than Ca. 



